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SUMMARY 


A generalized, real-time, piloted, visual simulation of a single-rotor 
helicopter, suspension system, and external load is described and validated for 
the full flight envelope of the U.S. Army CH-54 helicopter and cargo container 
as an example. The mathematical model described uses modified nonlinear clas- 
sical rotor theory for both the main rotor and tail rotor, nonlinear fuselage 
aerodynamics, an elastic suspension system, nonlinear load aerodynamics, and a 
load-ground contact model. The implementation of the mathematical model on a 
large digital computing system is described, and validation of the simulation 
is discussed. The mathematical model is validated by comparing measured flight 
data with simulated data, by comparing linearized system matrices, eigenvalues, 
and eigenvectors with manufacturers' data, and by the subjective comparison of 
handling characteristics by experienced pilots. A visual landing display sys- 
tem for use in simulation which generates the pilot's forward-looking real- 
world display is discussed, and a special head-up, down-looking load/landing 
zone display is described. 


INTRODUCTION 

Simulation and analytical research has been conducted at the Langley 
Research Center (LaRC) to compare various control system concepts for improving 
the handling qualities of single-rotor helicopters carrying relatively large 
external sling loads. These concepts include control jets at the load, a mov- 
able cable attachment point on the helicopter, and cable angle feedback into 
the helicopter stability augmentation system. It was believed that the most 
cost-effective and safe way to compare and study these systems was through the 
use of a piloted visual" simulation in which wide variations in parameters and 
concept optimization could be explored easily and quickly. 

Literature searches and discussions with industry did not locate any simu- 
lations for a helicopter sling load or any mathematical models having full- 
flight-envelope capability for a helicopter and load. Therefore such a model 
and simulation had to be developed. 

The purpose of this report is to describe the mathematical model of a heli- 
copter and external sling load which was developed, to describe the implementa- 
tion of the model to obtain a piloted visual simulation, and to describe the 
validation of the model and simulation through comparison of simulated data with 
actual flight data by means of analytical techniques and experienced pilots’ 
evaluations. 

The overall mathematical model is made up of numerous submodels that 
describe various components of the total dynamic system. The scope of these 
submodels is described in general terms as follows: 



Atmospheric model - The atmospheric model has variable air density, 
winds with variable magnitude and direction, and variable-intensity 
turbulence * 

Control system model - The helicopter control system model converts 
pilot* s cyclic-control-stick, collective-stick, and pedal inputs into 
main- and tail-rotor cyclic and collective-pitch inputs. 

Rotor models - Nonlinear models for the main and tail rotors define 
thrust, drag, and side forces as well as hub and force moments representa- 
tive of articulated rotors over a wide range of airspeeds from 100 knots 
through hover to rearward and sideward flight to at least 20 knots. The 
rotor models account for variable inflow velocity, variable rotor speed, 
blade twist, tip loss, blade coning, blade flapping, actuator dynamics, 
flapping-hinge offset, and tail-rotor flapping-hinge cant angle. 

Automatic flight control system (AFCS) model - The helicopter AFCS 
model is based on the system used in the U.S. Army CH-54 helicopter. This 
AFCS provides helicopter rate and attitude stabilization in roll, pitch, 
and yaw. 

Fuselage aerodynamics model - The fuselage aerodynamics model defines 
nonlinear lift, drag, and side forces as well as pitching, rolling, and 
yawing moments in terms of a wide range of fuselage angles of attack and 
sideslip, rotor down wash, body angular rates, and dynamic pressure. 

External-load aerodynamics model - An external-load aerodynamics 
model defines nonlinear lift, drag, and side force as well as pitching, 
rolling, and yawing moments as a function of angles of attack and side- 
slip, dynamic pressure, rotor downwash , and body angular rates. 

Load suspension model - A load suspension model defines cable tension 
in one or more cables and the resulting forces and moments acting on the 
helicopter and external load. 

Load-ground contact model - A load-ground contact model determines 
the approximate forces and moments acting on the load as it comes in con- 
tact with the ground for pickup and landing. 

A general set of nonlinear, rigid-body equations of motion for both the 
helicopter and external load determines the motion of each vehicle with respect 
to a flat, nonrotating Earth. An algorithm determines the trimmed helicopter 
control positions, helicopter attitude, and load position and attitude so that 
the entire dynamic system is in unaccelerated flight for a specified initial 
flight condition. Another algorithm obtains the equivalent linear system from 
the nonlinear model once the helicopter is trimmed; the linear system is used 
for verification and validation only. 

The mathematical model is programed on a general-purpose, real-time digital 
computer, with data for the U.S. Army CH-54 helicopter used as inputs; and 
appropriate outputs are fed to a cockpit having a set of flight instruments. 

The computer . outputs also drive a real-world, out- the -window visual display as 
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well as a computer-generated load/landing zone display for pilot use. The 
research pilot is able to control the simulated helicopter by making cyclic- 
and collective-stick and pedal inputs in the cockpit which generate electri- 
cal signals that are transmitted to the computer. Finally, provisions are made 
for recording simulated flight data and for interacting with the simulation 
from a control console. 

Use of trade names or names of manufacturers in this report does not con- 
stitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 

Contributions to this work and assistance to the authors were provided by 
the following persons: Dean E. Cooper, of Sikorsky Aircraft Division of United 

Technologies Corporation; NASA test pilot Perry L. Deal; Lt. Col, W, L. Welter, 
of Langley Directorate, USAAMRDL; pilots from Evergreen Helicopter, Incorpo- 
rated, McMinnville, Ore., Colonial Helicopters, Incorporated, Norfolk, Va. , 
and the 355th Aviation Company at Fort Eustis, Virginia; W. F. Lovell, J. B. 
Leavitt, and L. E. Becker, of Sperry Support Services; M. D. Par due, graduate 
student at Old Dominion University; and Lawrence E. Barker, Jr. , Lemuel E. 
Meetze, and Richard E. Bardusch, of LaRC. 


SYMBOLS 

Measurements, calculations, and programing were made in U.S. Customary 
Units. They are presented here in the International System of Units (SI) , 

Ajc/Bic rotor lateral and longitudinal cyclic control commands, 

defined by equations (12) , rad 

A ICaf cs ' B ICaf cs rotor lateral and longitudinal cyclic AFCS commands, 

defined by equations (14) and (13) , rad or deg 

A IC' B IC rotor lateral and longitudinal cyclic control displace- 

ments, defined by equations (23) , rad 

a rotor blade lift-curve slope, per rad 

a* small angle that defines the rotor drag force, defined by 

equation (35), rad 

a -i , a , a helicopter body-axes accelerations, m/sec^ 

x , n y , n z , n 

a Q rotor coning angle given by equation (29) , rad 

a l'kl rotor longitudinal and lateral flapping angles with respect to 

control axes, defined by equations (32) , rad 

aigrbfs rotor longitudinal and lateral flapping angles with respect to 

shaft axes, defined by equations (43), rad or deg 

B rotor-blade tip-loss constant 
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b 

C h' C £ 


number of blades per rotor 

Euler angle transformation matrix for helicopter and load, 
defined by equation (1) 

rotor torque coefficient, defined by equation (36) 
rotor thrust coefficient, defined by equation (27) 
rotor side-force coefficient, defined by equation (38) 
rotor blade chord, m 


>,d 


1c , e ' 2c , e ' 3c , e 


,a. 


cable direction cosines defined by equations (73) 


e kf' e kt 


mr 


fuselage and tail angle-of -attack corrections due to rotor 
downwash, rad 

rotor downwash factor, defined by equation (52) , rad 
rotor flapping-hinge offset, m 


^ F x,h' ^ F y,h' ^ F z,h' 
ZF x,£' EF y,£' ZF z,£ 
G B0' G Bq' G Bx' G A(j ) ' ] 

G Ap' G Ax' G 0tlj;' G etr' 

G e C h 

G gov ' K dgov ' K gov ' ^ 


force summations along helicopter and load body axes 
including all external forces due to rotor, body 
aerodynamics, ground contact, and suspension system, N 


AFCS feedback gains (see table I for units) 


G u' G v' G w 


g 

h h' h £ 

I bm' I bt 

^mr 


L pt 


engine/governor parameters (see table I for units) 
gust variables defined by equations (11) , m/sec^ 
acceleration of gravity, m/sec^ 
helicopter and load altitude, m 

mam- and tail-rotor blade flapping moment of inertia, kg-m 

2 

mam-rotor polar moment of inertia, kg-m 
engine-power- turbine moment of inertia, kg-m^ 


^xx^^yy,]!' ^zz,Yi r 1 helicopter and load mass moments of inertia about 

/ body axes, kg-m^ 

^•xx, £' -^yy , Z f ^zz, £ 
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I X2,h '^xz,# 

it 


helicopter and load products of inertia, kg-rr/ 

flow incidence at helicopter horizontal tail, defined by 
equation (54) , rad or deg 


ito fixed incidence of helicopter horizontal tail, rad 

K cO' K cl' K c2' K c3' 


K c4' K c5' K c6' K c7 


control-system stick gain, control mixing, and unit 
conversion constants (see table I for units) 


K 


fe 

k p 

K sc 

K v 


fuselage pitching-moment parameter due to rotor thrust used in 
equations (58) , N-m 

load-ground contact position force parameter used in 
equation (61), N/m 

load suspension system cable spring rate, N/m 

load-ground contact velocity force parameter used in equation (61) , 
N/ (m/sec) 


L c,£' M c,£' N c,£ 


L • 


ci , £ ' M ci , £ ' N ci , £ 


L d,h' M d,h' N d,h 


L d, £'%,£'%,£ 
L f,h' M f,h' N f,h 
^h'°h ,Y h 

L hub,h' M hub,h' N hub,h 
L hub , s ' M hub , s ' N hub , s 

L £' M £' N £ 

L £ ,D £ ,Y £ 


total load-ground contact moments, defined by 
equations (66) , N-m 

load moments due to ground contact at point i, N-m 

helicopter body-axes moments due to body angular rates, 
defined by equations (89) , N-m 

load moments due to load angular velocities, defined by 
equations (95), N-m 

total moments acting on helicopter due to fuselage aero- 
dynamics, defined by equations (58), N-m 

helicopter fuselage lift, drag, and side forces, N 

helicopter moments in body axes due to rotor moments 
transmitted at the hub, defined by equation (45) , N-m 

helicopter moments in shaft axes due to rotor moments 
transmitted at the hub, defined by equations (42) 
and (44) , N-m 

total aerodynamic moments acting on the load, defined 
by equations (96) , N-m 

lift, drag, and side forces acting on the load, defined 
by equations (92) , N 
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L r ^,M r ^ r N r ^ total body-axes moments acting on helicopter due to 

rotor forces and hub moments, defined by equation (46) , 
N-m 

L t h' M t h fN t h body-axes moments acting on helicopter due to suspension 

/ cable, defined by equations (75) , N-m 

L t £' M t £' N t £ body-axes moments acting on the load due to suspension 

cable, defined by equations (77) , N-m 

L u ,L v ,L w scales of turbulence, defined by equations (5) to (7) , m 

L wt ,M wt ,N wt body-axes aerodynamic moments acting on fuselage, 

appearing in equations (58) , N-m 

ZL^, E M h ,£N^ external moments, including main- and tail-rotor moments, 

fuselage aerodynamic moments, and cable suspension 
moments, N-m 

£ c suspension-system instantaneous cable length, defined by 

equation (69) , m 


suspension- system unstretched cable length, m 

mass moment of rotor blade appearing in equations (42) , kg-m 


m h' m & 1 

n l' n 2 < • • • ' n 6 


Ph^h' r h 'Vv*v r si 


helicopter and load mass, kg 


white-noise signals used in equations (10) 

helicopter and load angular velocity components in body 
axes, rad/sec 


P S 'V r 8 


helicopter angular velocity components in shaft 
axes, rad/sec 


P fQ t r 


helicopter angular velocity components in control 
axes, rad/ sec 


aerodynamic torque acting on main rotor, positive in direction 
opposite to rotation, N-m 

aerodynamic torque acting on tail rotor, positive in direction 
opposite to rotation, N-m 


shaft engine torque acting on rotor and fuselage (positive value 
tends to accelerate rotor and cause fuselage to yaw right ) , N-m 


Qg en gas-generator torque used in equations (48) , N-m 

Q s resulting main- or tail-rotor shaft torque acting on fuselage, N-m 
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dynamic pressure at fuselage and load, given by equations (50) 
and (91) , N/m 2 


R main- or tail -rotor radius , m 

S shaft tilt transformation matrix given by equation (18) 

s Laplace operator, s ec~^ 

T,H, J main- or tail-rotor thrust, drag, and side force defined by 

equations (28), (34), and (39), N 

T c suspension -system cable tension, defined by equation (70) , N 

u as h fV as h' w as h helicopter body-axes components of velocity, defined by 

equations (20) , m/sec 


u as , £ 9 v as , £ 7 w as , £ 
u cg , h ' v cg ,h ,w cg,h 
u cg , £' v cg , £ 9 w cg , £ 

u ci , e ' v ci ,e' w ci,e 


load body-axes components of total airspeed, defined by 
equations (90) , m/sec 

body-axes components of inertial velocity of helicopter 
c.g. , given by equations (80), m/sec 

body-axes components of inertial velocity of load 
c.g. , given by equations (80) with h replaced 
by £, m/sec 

Earth-axes components of inertial velocity of the ith 
load-ground contact point, defined by equation (60) , 
m/sec 


u hc g , e 9 V hcg , e ' w hcg , e 


u £cg , e 9 v £cg , e ' w £cg , e 


Earth-axes components of inertial velocity of 
helicopter c.g. , given by equation (84), m/sec 

Earth-axes components of inertial velocity of load c.g., 
given by equation (84) with h replaced by £, m/sec 


u gust , h f v gust , h 9 w gust , h ' 


u gus t , £ ' v gus t , £ 9 w gus t , £ 
u r,h' v r,h' w r,h 


helicopter and load body-axes components of gust 
velocity, given by equations (90) , m/sec 

body-axes components of total airspeed at rotor hub, 
given by equations (19) , m/sec 


u s' v s' w s 


shaft-axes components of total airspeed at rotor hub, 
given by equation (17) , m/sec 


Hw' V w /W w 


u wind , e ' v wind , e ' w wind , e 


control-axes components of total airspeed at rotor hub, 
given by equations (22), m/sec 

Earth-axes components of windspeed, given by 
equation (4) , m/sec 
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U wind , h ' V wind , h ' W wind , h ' l body-axes components of windspeed for helicopter 


U wind , £ ' V win d , £ ' W wind , £ 
V as,h' V as,£ 


and load, given by equation (4) , m/sec 


helicopter and load total airspeed defined by equa- 
tions (51) and (91) , m/sec 


V wind 


magnitude of wind vector, m/sec 




load sliding velocity beyond which sliding friction becomes inde- 
pendent of sliding velocity, m/sec 


w. 


as,£ 


load body z-axis component of total airspeed with constant added 
to improve numerical stability, m/sec 


X . , Y . ,Z . 

cx,e ci,e cx,e 


Earth-axes components of load-ground contact forces at 
the ith contact point, defined by equations (61) , 

(62) , and (63) , N 


X . n , Y . of Z . o 
cx,£' Ci,r cx,£ 


load body-axes components of load-ground contact forces 
at the ith contact point, given by equation (64) , N 


X c,£' Y c,£' Z c,£ 


load body-axes components of total load-ground contact 
forces, given by equations (66) , N 


x f,h' Y f,h' z f /h 


helicopter body-axes components of fuselage lift, drag, 
and side force, given by equation (57) , N 


X £' Y £' Z l 


lift, drag, and side forces expressed in load body axes 
(eq. (93)), N 


X r,h' Y r,h' Z r,h 


helicopter body-axes components of rotor forces, given 
by equation (41), N 


X r,s' Y r,s' Z r,s 


shaft-axes components of rotor forces, defined by 
equations (40) , N 


X t,h' Y t,h ,Z t,h 


helicopter body-axes components of cable tension, given 
by equation (74) , N 


X t,£ ,Y t,£' Z t,£ 

x a,h' y a,h' z a,h 


load body-axes components of cable tension, given by 
equation (76) , N 

helicopter body-axes distances from helicopter c.g. to 
cable attachment point, m 


c a,£' ^a,£ ,z a,£ 


load body-axes distances from load c,g. to cable attach- 
ment point, m 


x ci , e ' ^ci , e ' z ci , e 


Earth-axes coordinates of ith ground contact point, m 
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x ci,£' y ci, V z ci, £ 
X col ' x lat 9 x lon 9 X ped 


< ha / e ' y ha , e ' z ha, a 


ncg ,e r y hcg , e 9 z hcg , e 


c £a , e 9 y £a , e ' z £a , e 


c £cg , e ' y £cg , e 9 z £cg , e 


C £cg , h 9 ^£cg , h' Z £cg , h 


X £p , h 9 ^£p ,h' Z £p,h 


x £p,£' y £p, £ #z £p, £ 

x mr , h ' ^mr , h ' z mr , h 9 
X tr,h/ y tr,h' z tr,h 
x ps , h ' y ps , h ' z ps , h 


x screen ' y screen 


x v,e' y v,e' z v,e 
x v,h ,y v,h' z v,h 
X wt,h'^wt,h' Z wt,h 
x zp , e ' y zp , e ' z zp , e 


load body-axes distances from load c.g. to ith ground 
contact point, m 

control displacements of collective stick, lateral and 
longitudinal cyclic stick, and pedals (positive dis- 
placements cause climb, right roll, pitch-up, and yaw 
left, respectively), m 

Earth-axes inertial position of helicopter cable attach- 
ment point, given by equation (67) , m 

Earth-axes inertial position of helicopter c.g. , given 
by equations (85) , m 

Earth-axes inertial position of load cable attachment 
point, given by equation (68) , m 

Earth-axes inertial position of load c.g., given by 
equations (85) with subscript h replaced by £, m 

helicopter body-axes coordinates of load c.g., given by 
equation (97), m 


helicopter body-axes coordinates of any point p on the 
load, given by equation (98) , m 


load body-axes coordinates of any point p on the 
load , m 


helicopter body-axes coordinates of main- and tail-rotor 
hub, m 

helicopter body-axes coordinates of eye-level position 
of pilot in left seat, m 

nondimen sional screen coordinates of any point p, 
used in load/landing zone display defined by equa- 
tions (101) , m 

Earth-axes coordinates of helicopter viewpoint, given 
by equation (99) , m 

helicopter body-axes coordinates of down-looking 
viewpoint, m 

helicopter body-axes coordinates of wind-tunnel mounting 
point, m 

Earth-axes inertial position of any point p fixed on 
the Earth, m 
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x zp,h' Y zp,h' 


a , 3 
c , e c , e 


a f l 

a £' g Jl 

3 


z zp 7 h helicopter body-axes coordinates of any point p on 

' the Earth, m 

cable angles with respect to the vertical, defined by equations (71) 
and (72) , rad 

fuselage angle of attack and sideslip, defined by equations (49) 
and (55) , rad or deg 

local fuselage angle of attack, defined by equation (53) , rad or deg 
load angles of attack and sideslip, defined by equations (91) , rad 
rotor orientation angle, rad 


Y 


rotor Lock number. 


pacR 


Rotor blade flapping moment of inertia 
ALfh, ^ L 2h incremental fuselage lift, N 

incremental fuselage pitching moment, N-n 


-m 


At 

6 


3t 


0 


ct 


7 0m 


Ot 


simulation integration step size, sec 

tail-rotor flapping-hinge cant angle, rad 

commanded value of tail-rotor collective pitch given by 
equation s ( 12 ) , rad 

main-rotor collective pitch, given by equations (12) , rad or deg 

effective value of tail-rotor collective pitch, given by 
equation (47) , rad or deg 


0 f main- and tail-rotor AFCS collective pitch command, defined 

by equations (15) and (16), rad 

0 S longitudinal shaft tilt angle, positive for thrust vector tilted 

from — z^ toward -x^, rad 


0 


.75 


h'k'k 


rotor blade twist angle from root to tip, rad 

rotor collective pitch at three-fourths radius, defined by 
equation (33) , rad 

helicopter pitch, roll, and yaw AFCS attitude errors, rad 
main- and tail-rotor inflow ratio, defined by equation (25) 
coefficient of sliding friction used in load-ground contact model 
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V 

p 

Pa 

o 


T eng 

x 63t 


t A 

WV 


main- and tail-rotor tip-speed ratio, defined by equation (24) 
rotor induced inflow ratio, defined by equation (26) 
atmospheric density given by equation (2) , kg/m 
swashplate actuator damping ratio 

rotor solidity, — 

TTR 

body-axes components of rms gust intensity, m/sec 

engine time constant, sec 

tail-rotor time constant, sec 

rotor inflow time constant, sec 

white noise power spectral density, sec 

roll .(<j>) , pitch (0) , and yaw Op) Euler angles of 
helicopter and load, rad 


*3 


^wind 


lateral shaft tilt angle, positive for thrust vector tilted from 
-z^ toward +y^, rad 

wind angle measured clockwise from true north that defines the 
direction from which the wind is coming, rad 


wind-tunnel yaw angle, rad 
ft rotor angular velocity, rad/sec 

ft D commanded rotor angular velocity, rad/ sec 

ftp t power-turbine angular velocity, rad/sec 

03 a swashplate-actuator natural frequency, rad/sec 

Abbreviations : 

ADC analog to digital converter 

AFCS automatic flight control system 

c.g. center of gravity 

DAC digital to analog converter 
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BLCG 

buttock line at c.g. 

FSCG 

fuselage station at c.g. 

WLCG 

water line at c.g. 

RTS 

real-time system 

rms 

root mean squared 

VLDS 

visual landing display system 


Dots over a symbol denote derivatives with respect to time. 

Initial values of variables are denoted by (G) following the variable 
symbol. 

Matrix transpose is denoted by a superscript T. 


DESCRIPTION OF MATHEMATICAL MODEL 

The mathematical model for the simulation of a single-rotor helicopter and 
external load can be given in terms of submodels for the various components of 
the total dynamic system. 

A block diagram for the overall mathematical model is given in figure 1. 
This diagram together with the Contents should be useful in understanding the 
interrelationship between the submodels. In the following sections the mathe- 
matical model is described in approximately the order it would be executed in 
the computer so that the reader can write his own simulation program more 
easily. 


Coordinate Systems 

The following right-hand orthogonal axis systems are used in this report: 

(1) Earth axes (subscript e) : origin fixed on the Earth* s surface, 

x e -axis pointing north, y e -axis pointing east 

(2) Helicopter body axes (subscript h) : origin at the helicopter c.g. , 

x^-axis pointing forward in the plane of symmetry of the fuselage 
and parallel to the helicopter water line, z^-axis pointing downward 
away from the main rotor and in the plane of symmetry (see figs. 2 
and 3) 

(3) Shaft axes (subscript s) : origin at the center of the rotor hub, 

y s -axis rotated through the lateral shaft tilt angle <j) s about the 
x^-axis, x s -axis rotated through the longitudinal shaft tilt 
angle 0 g about the y s -axis, z s -axis coincident with the rotor 
shaft, applies to both the main rotor and the tail rotor (see fig. 3) 
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(4) Control axes (subscript w) : origin at the center of the rotor hub, 

z w -axis directed toward the fuselage along the axis of no feather- 
ing (the physical axis of a pure flapping rotor - a rotor with 
blades fixed in pitch but free to flap) , x w -axis chosen to give no 
y w -component of velocity relative to the free-stream air (the free- 
stream air to include the air motion from turbulence and steady 
winds) r applies to both rotors (see fig. 3) 

(5) Wind-tunnel axes (subscript wt) : origin at the helicopter or 

external-load wind-tunnel mounting point, x wt -axis pointing into 
the relative wind, z wt -axis pointing downward and perpendicular to 
the wind 

(6) External- load body axes (subscript £) : origin at the load c.g. , 

x^-axis pointing forward in the plane of symmetry of the load, 
z^-axis pointing downward in the plane of symmetry of the load 

The Earth- to-body axis transformation for the helicopter based on the 
standard yaw, pitch, and roll Euler angle rotation sequence shown in figure 3 
is given in matrix form as 


-h = 


cos 9^ COS 
sin 0^ sin (f>^ cos 
- cos (J)^ sin 
sin 0^ cos (j)^ cos ij)^ 
■+ sin cj)^ sin ^ 


cos 

0^ sin 



-sin 

A 

CD 


sin 

CD 

sin <£ h 

sin 


cos 0^ 

sin 


+ 

cos 

4> h cos 





sin 

e h 

COS (f)^ 

sin 


cos 0^ 

cos 



- sin <1).^ cos ^ 


CD 


The transformation matrix for the load has the same form with the subscript h 
replaced by subscript £. 


Atmospheric Model 

The atmospheric model used allows for air-density variations with alti- 
tude, variable wind direction and magnitude, and variable-intensity atmospheric 
turbulence. The air density is calculated according to the polynomial in h 


p = 1.2266 - (1.176 x 10" 4 )h + (4.337 x 10 9 )h 2 - (7.463 x 10“ 14 )h 3 

+ (5.538 x 10“ 19 )h 4 - (9.357 x 10~ 25 )h 5 (2) 

where h is the helicopter or external- load altitude above sea level, in 
meters, as determined from integrating the equations of motions discussed below. 
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The wind 


Steady winds are specified in terms of magnitude and direction . 
is expressed in Earth axes by the relations 

u wind,e ~ V wind cos ^wind 
v wind,e “ V wind s ^ n ^wind / 
w wind , e = 0 


(3) 


where ^ w i nc j is ^e angle measured clockwise from true north that defines the 
direction from which the wind is blowing. Steady vertical winds are not simu- 
lated. The components of the steady winds are expressed in helicopter body 
axes by using the transformation 


r > 

r \ 

^ind , h 


^indye 

v wind , h 

) - c h< 

v wind,e 

w wind,h 


w wind,e 

V. J 

v J 


(4) 


The components of the steady winds are expressed in load body axes by using 
equation (4) with the subscript h replaced by subscript £. The atmospheric 
turbulence mathematical model is based on the Dryden spectrum of turbulence as 
discussed in reference 1. This theory allows the simulation of atmospheric 
turbulence by passing uncorrelated white noise through linear filters to obtain 
gust components in body axes. 


The scales of turbulence for the helicopter and the load are calculated 
according to the relations shown for the helicopter 


65.7h h 

L u,h = { 533 - 4 
762.0 


1/3 


v 


,h ~ L u,h 


r 


h v 


L w,h = ( 533 - 4 


762.0 




(h h ^ 762.0; O n < 6 .4) 
(O n > 6 . 4 ) 
(h h > 762.0; O u < 6 .4) 




(h h £ 762.0; O u < 6.4) 

(a u ^ 6.4) ) 
(h h > 762.0; a u < 6.4) 


(5) 


(6) 


(7) 
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where a u is the longitudinal rms gust intensity (assumed to be the same for 
the helicopter and load). The scales of turbulence for the load are calculated 
by using equations (5), (6), and (7) with the subscript h replaced by sub- 
script A, The vertical and lateral gust intensities for the helicopter are 
given by 


J w,h 


0.0376a u h h 


1/3 


and 


a v = a u 


(a u a 6.4) 
(a u < 6.4) 


The gust intensities for the load are calculated from equations (8) with the 
subscript h replaced by subscript A. The white noise power spectral density 
is given by 

= At/2TT 


where At is the digital simulation computation step size in seconds. 

The components of gust velocity for the helicopter and load in body 
axes are calculated by passing six uncorrelated white noise signals 
n l' n 2 f • • •/ n 6 through linear filters and are given as follows for the 
helicopter: 


u gust,h 


"u,h 


s + 


v as,h 






G v,h 


s + 


V 


gust ,h 


s + 


v as,h 

jf£ L v,t 

v as , h \ 
L v,h j 


n 2 


( 10 ) 


J w,h 


s + 


w 


gust ,h 


s + 


V as,h 
h Wf h 
v as f h \ 
L w,h / 
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where 


G u,h = a J V as,h /(^o L u,h) 

G v,h = °v \J^ V as ,h /( 2Tr $ o L v,h^ ^ 

G w,h = v 


(ID 


The gust components for the load are computed in a similar manner. The filters 
may be implemented digitally by using Z transform techniques. 

It was judged after discussion with G. H. Fichtl, of the NASA Marshall 
Space Flight Center , that for sling-load cable lengths of approximately 
30 meters, the gusts for the helicopter and the external load will be essen- 
tially uncorrelated. Thus two sets of uncorrelated white noise signals are 
used for the helicopter and the external load. 


Control System 


The flight control system mathematical model converts the pilot* s cyclic- 
stick, pedal, and collective- stick displacements and AFCS outputs to equivalent 
main-rotor and tail-rotor control inputs. The main-rotor collective pitch 0 Qm , 
longitudinal and lateral cyclic pitch B IC and A IC , and tail-rotor collective 
pitch 0 ct are given in terms of collective, cyclic, and pedal displacements 
and AFCS outputs by the following expressions obtained from Sikorsky Aircraft 
Division of United Technologies Corporation: 


®0m K c0 + K cl x col + ®mafcs 
B IC = K c2 x lon + B ICafcs 

> ( 12 ) 

A IC = K c3 x col + K c4 x lat + A ICafcs 
®ct = k c5 + K c6 x ped + K c7 x col + ®tafcs 


where x col , x lon , x lat , and Xp ed are the pilot's control displacements 
with respect to a fixed reference position. The constants K c0 through K c 7 
are used for stick gain, control mixing, and unit conversions. Since control 
feel is important to pilots, a programmable hydraulic control loader is used 
for the cyclic control stick and the pedals. Pertinent details of this system 
are discussed in the section "Simulation Description." 
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Automatic Flight Control System Model 


The mathematical model of the automatic flight control system (AFCS) is 
similar to that of the attitude command system employed in the U.S. Army CH-54A 
helicopter, and was obtained from unpublished Sikorsky Aircraft data. The 
pitch channel equation is given as 


B ICafcs ~ G B6®h + G Bq% + G Bx x lon 


(13) 


where 

e h - e h - e h (o) 

and 

x lon = x lon " x lon (0) 

The quantities 0^(0) and x i on (°) are the precalculated initial trim values 
of 0^ and x ]_ on discussed in the section "Trim Calculations.” 

The roll equation is given similarly as 


A ICafcs G A(j)^h " G Ap p h + G Ax X lat 


(14) 


where 


and 


= - ^h (0) 


x lat = x lat " x lat (0) 


The yaw channel is modeled by the following equation: 


A 

®tafcs “ G 0t^h + G 0tr r h 


(15) 


where 


A 
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The heading command i|^ q is initially set to the helicopter trim yaw angle . 
During the simulation whenever the pilot actuates either microswitch located on 
the simulator pedals, then ifj Q is set to the current value of When the 

pilot releases the micro switches, ip Q is set to the value of ijj^ just prior 
to switch release. If the heading hold feature is not desired, then a push 
button discrete in the cockpit will set the gain Gg^ to zero. 

The altitude channel is modeled by using the following equation: 


A. 

®mafcs ~ G 8ch^h 


(16) 


where 


A 



= h h " 


h 


o 


The altitude command h Q is initially set to the helicopter altitude. During 
the operate mode, the pilot can reset h Q to the current altitude or disable 
the altitude hold feature by using a push-button discrete in the cockpit. 


Rotor System 


In the interest of real-time pilot-in-the-loop computer simulation, every 
effort is made to keep the rotor mathematical model simple but adequate to 
allow forward flight to at least 100 knots, rearward and sideward flight to 
20 knots, autorotations, and large-angle maneuvers. The forces and moments due 
to the main and tail rotors are calculated by using a modified Bailey represen- 
tation discussed in references 2 to 5. 

The following discussion of rotor modeling is general and is applicable to 
both the main and tail rotors. Important limitations are given when necessary. 
Rotor variables will need "main (m) " or "tail (t) " identification when included 
in a computer program. 

The velocity of the free stream at the hub is expressed in shaft axes in 
terms of inertial velocity of the helicopter center of gravity, body angular 
rates, steady winds, atmospheric gusts, and position of the main- and tail- 
rotor hubs by the following expression: 


c 

r 's 

u 


u , 

s 

L i 

r,h 

V 

\ = S / 

v , 

l s 

/ \ 

r,h 

w s 


w r,h 

V J 


v J 


( 17 ) 
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where 


S = 


cos 0 Q 
0 


" sin 9 S 


cos (f> s 


cos 0 S 


sin (f) s 


-sin <P s 
cos 0 S cos $ s 


(18) 


and 


u r,h = u as,h + % z r,h 


v , - v , + r, x -i 

r f h as,h h r,h 


w. 


r,h w as,h + ^h^r,h 




r h^r ,h 
p h z r,h ) 
q h x r f h 


(19) 


and 




U , = U , + U . , + U . J, , 

as,h cg,h gust, h wmd r h 


v as,h v cg,h + v gust,h + v wind,h 7 


w as,h w cg,h + w gust,h + w wind,h 


( 20 ) 


The rotor orientation angle 3 is defined by the relation 


3 = tan' 


-1 v s 


(21) 


and the velocity of the free stream at the hub is expressed in control axes by 
the approximation given in reference 2 (see fig. 3(c)): 




= U s cos 3 + v s sin 3 


v w = 0 


w w - w s - B IC u s " A IC v s 


( 22 ) 
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where B^ c and A-J- c are the shaft relative cyclic control inputs for the main 
rotor * For the tail rotor , A® and B* are of course zero. 


Actuator dynamics are modeled by passing the control inputs through second- 
order filters as 


IC " S 2 + 2p a U) a s + lo a 2 IC 


IC s 2 + 2p a i0 a s + » a 2 IC 


(23) 


The filters may be implemented by using Z transform techniques. The cyclic 
control inputs A IC and B IC are generated by pilot cyclic stick motions and 
automatic control system output given by equations (12) to (14). 

The rotor forces and moments are functions of the rotor tip-speed ratio 
and induced inflow ratio, as discussed in references 2, 4, 6, and 7. The tip- 
speed ratio y is given by 


U = 


^w 

ftR 


(24) 


where ft is the rotor angular velocity. The inflow ratio A is calculated by 
the implicit equation 


A = 




w 




- V 


(25) 


where the induced inflow ratio V is determined from the differential equation 



This method of computing the inflow ratio A assumes that the inflow is con- 
stant across the rotor disk. The thrust coefficient C T is calculated below, 
and a value of the time constant of approximately 0.1 sec is chosen to 

approximate the time lag associated with change in rotor inflow. This tech- 
nique of using a first-order differential equation to calculate the induced 
velocity ratio V is superior to algebraic calculation, because algebraic 
methods were found to be numerically unstable in digital simulations. 
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The rotor thrust T and coning angle ag are calculated to the third 
power of the tip-speed ratio according to the following relations taken from 
reference 4, where higher order terms have been neglected; 


C«p 


~ B 2 + ~ y 2 )A + 





+ 






( 27 ) 


where 0 q is the effective blade pitch angle at the blade root and 0]_ is the 
twist of the blade. Then 


T = bcRp(Rft) 2 




(28) 


and 


= Y 


B 3 + 0 . 04 y 3 


A + 


1 4 


1 7 2 

3 B y 


K + fe 


12 


B 3 y 2 


! ) 0 1 


(29) 


The term in equation (22) in reference 4 involving the blade mass moment con* 
tributes less than 0.5°, is essentially constant, and is neglected here. 

The fuselage angular velocity expressed in control axes is required in the 
calculation of certain rotor forces and moments. This quantity is obtained by 
rotating the fuselage angular velocity expressed in shaft axes through the 
rotor orientation angle 3, neglecting the small cyclic pitch angles A IC 
and B IC . (See fig. 3(c) for clarification.) 


p w = p s cos ^ + q s sin 0 
q w = _p s sin P + cos & 




J 


where 


(30) 


fp A 

F s 



q 5 

> = s T <: 

% 

< x sj 




(31) 
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The flapping angles a^ and are calculated relative to control axes 

by using the following formulas derived from references 3, 6, and 8 and unpub- 
lished data obtained from Sikorsky Aircraft Division of United Technologies 
Corporation: 


ai = 


1 - 


( 2 A + 3 © , 75^ y + 


P w 


8 


B 4 yfi 


2B^ 


and 


*>1 


,2 \3 Pa ° 


1 + 


% 16 *V 

& B 4 y8 j 


2B^ 


J 


(32) 


For a rotor blade with linear twist and constant chord it can be shown that 
replacing 9q, the blade pitch at the root (appearing in the references), with 
®.75' the pitch at three-fourths radius, and dropping 9^ will have a negli- 
gible effect on the overall solution. An expression for 0 ^ is given as 

9.75 * 9 o + °- 759 i (33) 

The downwind horizontal component of the rotor force in control axes is 
expressed as 


H = Ta’ 


(34) 


where the small angle a* is a function of the useful and induced rotor drag- 
lift power and inflow; however, it behaves similarly to the longitudinal 
flapping angle a^. An expression for a* which includes the influence due 
to body rate, as discussed in reference 8 and in unpublished Sikorsky Aircraft 
data, is as follows; 


■ 

a 



( 2A + 3 


^2»(l 

B 4 y8 \ 


°- 2 ".7 5 \ 
c T /a j 


(35) 


The expression for the rotor torque which accounts for both acceleration 
and deceleration torques is derived from equations (9) and (11) of reference 4 
and unpublished data obtained from Sikorsky Aircraft. The torque coefficient 
may be expressed as a polynomial in y times the major variables as 
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^2- = (o. 00109 - 0 . 00 36 A - 0.00270 ?5 - 1.10A 2 - O.545A0 ?5 + O.1220 2 ?5 ) 

+ (0.00109 - 0.00270 75 - 3.13A 2 - 6.35A0 ?5 - 1.930 2 ?5 )p 2 
- O.133A0 75 y 3 + (-0.976A 2 - 6.38A0 ?5 - 5.260 2 ?5 )y 4 (36) 

Thus, 


Q r = bcK p (Rfi) 


2 C Q 


(37) 


The torque acting on the main rotor Q is calculated by using main-rotor 
parameters in equations (36) and (37) . The reaction torque on the fuselage, 
which is a function of Q ., is calculated in the engine dynamics and governor 
model discussed in the next section. The torque acting on the tail rotor .Q ^ 
is calculated by using tail-rotor variables in equations (36) and (37) . It is 
assumed that the tail-rotor reaction torque acts on the fuselage directly and 
is equal to Q at * 

The rotor side force J in control axes is calculated from the following 
expression derived from equation (3) in reference 3, assuming that terms 
involving pitch and roll rate may be neglected and that the blade pitch 0 Q 
can be replaced by 0 ^ and terms involving 0j_ dropped: 


C Y 

a 


a_ 

2 


T b l X ~ 2 a 0^ X + 7 a l b l^ “ a 0 a l^ + t a 0 £ 


- (! i“o - T b i - I ' j2b i) e . 


75 


(38) 


from which 


J = bcRp(Rfi) 



(39) 


The rotor forces in control axes are resolved into shaft axes by assuming 
that the components of thrust along the x s - and y g -axes (fig. 3(b)) are TBj C 
and TA j£/ respectively, and that the components of drag H and side force J 
along the z s -axis are negligible compared with the thrust, so that 


X r s = ~H cos 0 - J sin 0 + TB ^ 

Y r s = " H sin 2 + J cos 0 + TAj C > 


(40) 


Z 


r,s 


-T 


J 
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These forces are expressed in helicopter body axes by using the transformation 


fX_ V, 

r,h 


Cy. ~\ 
x r,s 

Y . 

\ = s / 

Y 

r ,h 

s 

r,s i 

< Z r,hj 


< z x , s J 


(41) 


The hub moments due to flapping hinge offsets are given in shaft axes by 
the following approximate expressions derived from reference 3, where higher 
order terms have been neglected: 


L hub , s 2 I "V D ls 


Mi 


hub, s 


ebft M^c 


Is 


(42) 


where the cyclic flapping angles relative to shaft axes are given by 


a = a 1 cos 3 + sin 3 - B ^ 




IC 


and 


bfs = b-^ cos 3 ~ sin 3 + 


IC 


(43) 


The torque about the rotor shaft is given by 


N hub,s ®s 


(44) 


For the main rotor, Q s is equal to Q eng , the engine torque computed in the 
section "Engine Dynamics and Governor Model." For the tail rotor, Q s is 
equal to These hub moments are expressed in helicopter, body axes by the 

trans forma txon 


C L > 

hub , h 


Xiub , s''' 

\ M hub,h 

> = S < 

M hub,s , 

V N hub,h> 


^hub , s> 


(45) 


24 



The total moments acting on the fuselage due to the main rotor and tail 
rotor are calculated by using the rotor forces given by equation (41) and hub 
moments given by equation (45) as 


r T 

r,h 


L hub,h 


0 

i 

V 

y r,h 


^ X r,h^ 

M r,h 

>■ < 

M hub , h 


*r,h 

0 

_x r,h 

< 

Y r,h 

i 

v N r,h> 


l N hub,h^ 


Z*r, h 

x r ,h 

0 


^ Z r,h> 


(46) 


The mathematical model developed here assumes no 63 ^ hinges on the main 
rotor, that is, blade-coning and lead-lag motion does not affect blade pitch. 

It is assumed, however, that the tail rotor has large values of 63 ^; for 
example, if 6 3t = 45 °, then 1 ° of increased coning reduces the blade pitch by 
1° and conversely. This effect is modeled as follows: the tail-rotor coning 

angle ag t is calculated using equation (29) with the current value of effec- 
tive tail-rotor collective pitch 0Qt* The new value of 0 q^ is determined by 
solving the first-order differential equation 


5 0 t 


T <53t (9ct 


- a Qt tan $ 3t 


V 


(47) 


where 0 Q t is the collective pitch value commanded by the pilot's pedals and 
the automatic control system. The time constant is taken as small as 

possible while maintaining a good margin of numerical stability. A value 
between 0.05 and 0.2 second is typical. This method is superior to a purely 
algebraic calculation of 0 Qt in that the algebraic method was found to be 
numerically unstable in digital simulations. 


Engine Dynamics and Governor Model 

The engine dynamics and governor model was adapted from one used by Boeing 
Vertol Company in a simulation of a heavy- lift helicopter. This model includes 
the effects of a gas generator, a power turbine, rotor inertia, and shaft com- 
pliance. The following differential equations define this model: 

^m “ [peng ~ ^am + K dgov^pt 
Seng = ~ 

V * [Bgen + K gov< n o " V " 

2 gen ” (pam “ 2 gen + G gov^o 


•V 1 




q ) /i 

m / mr 


Seng “ K dgov^pt ~ ^mji] pt 

X, 


> (48) 


^pt 5 


eng 
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The input variables to these differential equations are the reference rotor 

speed, and Q am , the aerodynamic torque acting on the main rotor, given by equa- 
tion (37) with the appropriate parameters for the main rotor. The outputs are 
the rotor speed ^ and the engine torque Q e . The power-turbine speed 
and the gas-generator torque Q gen are internal variables. The constant 
represents the shaft compliance. The constant ^ gov is needed for computa- 
tional stability. The constants K gov and G gov ^ are the power-turbine and 
gas-generator speed governor gains, and T eng is the gas-generator time con- 
stant. The aerodynamic torque Q am appears in the Q gen differential equa- 
tion to allow the governor to hold reasonably constant rotor speed under widely 
varying aerodynamic torques. 


Fuselage Aerodynamics 

The helicopter fuselage force and moment data are assumed to be given in 
equation and table form in wind-tunnel axes in terms of local angle of attack, 
local angle of incidence at the tail, angle of sideslip, body angular rates, 
and dynamic pressure. 

The free-stream angle of attack is given by 

-1 w as,h / < < v 

oif = tan x — — — ( — Tf .== a f = tt) (49) 

u as,h 

where the velocity components are calculated by using equations (20) . The 
free-stream dynamic pressure is given by 

5h - 2 P V L,h < 50 ' 


where 


V 

as 




i 


2 2 2 

u , + v , + w , 
as,h as,h as,h 


(51) 


The effect of the main-rotor downwash on the local angle of attack is 
accounted for by introducing a rotor downwash factor obtained from unpublished 
Sikorsky Aircraft data as 


e mr 


C 


Tm 




( 52 ) 


where A m and y m are the main-rotor inflow and tip-speed ratios, respec- 
tively, and Crj^ is the main-rotor thrust coefficient. The local angle of 
attack for the fuselage is given by 
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( 53 ) 


a f SL a f e mr e kf 


(-tt = a f £ = 7T) 


and the local incidence at the tail is given by 


i t “ 1 to " (e kt “ e kf )e mr 


( 54 ) 


where e^ t and e, ^ are empirical constants which have been determined by 
Sikorsky Aircraft from wind-tunnel and flight-test correlation. The constant 
i to is the fixed incidence of the horizontal- tail surface. 

The fuselage sideslip angle is given by 

3 f = sin {- - = 3f = 2 ) (55) 

as ,h 


and the wind-tunnel yaw angle is given by 


° -0f 


( 56 ) 


Since wind-tunnel data generally do not cover the full ranges of angle of 
attack and sideslip, it is assumed that force and moment coefficients remain 
constant beyond the limits of these angles. This assumption is based on the 
fact that generally when these angles are large, the airspeed is low, so that 
the fuselage forces and moments are relatively small. 

The forces in wind-tunnel stability axes are transformed into body axes by 
the relation 
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( 57 ) 


The basic fuselage aerodynamic moments are assumed to be given in body axes, 
and the total fuselage aerodynamic moments, including the effects of the c.g. 
being offset from the wind-tunnel mounting point as well as damping due to 
angular velocity and rotor downwash , are obtained as follows: 
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L f,h ~ H/t + Z f ,h Y wt,h “ Y f ,h z wt,h + 

M f,h = M wt + X f,h z wt,h " Z f,h x wt,h + M d,h + K fe T m ' ) 

N f,h _ N wt + Y f,h x wt,h ” X f,h Y wt,h + N d,h j 


( 58 ) 


The last term in the Mf ^ equation represents the static moment due to main- 
rotor downwash at the horizontal tail. The terms and h 

account for the aerodynamic moments due to body angular velocities and are dis- 
cussed subsequently in "Application to U.S. Army CH-54 Helicopter and Cargo 
Container, " 


External-Load Aerodynamics 

The aerodynamic forces and moments acting on the external load are calcu- 
lated in essentially the same way as the fuselage aerodynamics: in terms of 

dynamic pressure, angles of attack and sideslip, and body angular rates. The 
rotor downwash effects are not considered important for long cables as far as 
pendulum and rocking motions are concerned, but experience has shown that some 
simulated downwash is necessary to provide aerodynamic damping of vertical 
bounce-type oscillations near hover due to cable elasticity. Generally force 
and moment data are missing at large angles of attack, so that some type of 
trigonometric formulation for forces and moments is required. Specific details 
of the aerodynamic characteristics for a particular external load are discussed 
subsequently in the application section. 


Load-Ground Contact Model 


Since the pickup and release of sling loads are to be studied, a mathe- 
matical model of the load-ground contact forces is required. This model is 
derived on the assumptions that (1) the ground can be represented by springs 
and dashpots , so that vertical forces act on each corner of the base of the 
load in proportion to the distance that the corner is below ground level, and 
(2) viscous sliding friction retards translational motion. 

The inertial position of the ith contact point is determined from the 
relations 


f x ""n 

ci , e 
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Y Icq , e 
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v z 5/cg,e ^ 


^ z ci,£ > 


( 59 ) 
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where the first term on the right-hand side defines the inertial position of 
the load c.g„; the second term, the contact point with respect to the c.g* The 
inertial velocities of this load contact point are given by 


Cm . ^ 

ci,e 

v ci,e 
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u £cg,e 

V £cg,e 

> + c i t < 

q £ z ci, £ ~ r £ y ci,£ 
r £ x ci,£ - p £ z ci,£ 

^ w ci , e j 
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\ (i = 1, 2, 3, 4) 


(60) 


The vertical ground contact forces are computed for each corner of the load by 
the equation 


(i = 1, 2, 3, 4) (61) 


If Z ci e > 0, then the load is being pulled toward the Earth; therefore, 
Z c± e Is set to zero* The viscous sliding forces are then determined in 
terms of the vertical force and the horizontal component of velocity as 
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(63) 


< v ci f e < - v o> 


for i = 1, 2, 3, 4. This model assumes that the sliding force is propor- 
tional to the vertical force Z . for sliding velocities greater than 

Cl |6 

V Q = 0*31 m/sec and proportional to the product of vertical force and 
sliding velocity for velocities less than this value. The load contact 
forces at the ith contact point are transformed into load body axes by the 
transformation 
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(i = 1, 2, 3, 4) (64) 


and the moments acting on the load due to ground contact forces at the ith 
contact point are 


L ci, Z 
M ci , Z 
N ci,l 


Z ci,Z Y ci,Z 

x ci,Z z ci,Z 
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Y ci,Z z ci,Z 
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(i = 1, 2, 3, 4) 


The total forces and moments acting on the load are thus given by 


(65) 


X c,£ “ Zx ci,£ 

Y c,& “ EY ci,£ 

Z c,£ = Ez ci,Jl 

and \ (i = 1, 2, 3, 4) (66) 

L c,£ ~ ^ L ci,& 

%,i - 

a c,8, ■ Zl) ci,«. 

J 

Load Suspension System 

The mathematical model for the external-load suspension system is based 
on the assumption that the cable (s) may be represented by spring(s) without 
damping. This assumption allows the cable tension to be calculated easily in 
terms of the cable spring constant and the vector distance between the cable 
attachment points. With a digital simulation the inertial positions of the 
cable attachment points are readily calculated, thus simplifying the calcula- 
tion of the distance between the points. It was found by trial and error that 
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the simulated cable spring constant should be selected such that the cable 
vertical bounce frequency is not higher than approximately 2 Hz, in order to 
maintain numerical stability with an integration step size of 1/32 second in 
digital simulations. 

The inertial position of the cable attachment point on the helicopter 
shown in figure 4 is calculated in terms of the inertial position of the heli- 
copter c.g. and the distance between this position and the attachment point as 
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The inertial position of the cable attachment point on the load shown in fig- 
ure 4 is calculated similarly as 
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The cable length is determined from the square root of the sum of the squares 
of differences in inertial coordinates of cable attachment points as 


*c ]j( x ha r e “ x £a , e ) + (^ha f e ~ y £a,e) + ( Z ha,e ~ z Za, 


e) 2 < 69 > 


If the unstretched cable length is & c q, then the cable tension is simply 


T c K sc ^c ~ ^c(P 


(70) 


where K gc is the cable spring constant. A test on the sign of T c is con- 
tinuously made to insure that the cable does not exert a compression force. 

If (£ c - £ c0 ) becomes negative, then is set to zero. 

The cable tension force must be expressed in helicopter and load body axes 
so that its effect can be included in the equations of motion. Before the 
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tension force can be expressed in body axes, the inertial orientation of the 
cable - assumed straight - must be determined. 

The angle that the cable makes with the vertical in the north-vertical 
plane is defined in figure 4 and is given by 


8 . = tan" 1 

C ' e z £a,e 
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(71) 


The inertial cable angle in the east-west direction is also defined in figure 4 
and is given by 
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Three cable direction cosines are defined with respect to Earth axes by 

a 

d lc,e = sin 3 c,e cos a c,e * ( Ha,e ~ x ha,e>Ac 

d 2c,e = “ sin a c,e = ^la r e ~ / (73) 

d 3c,e = cos e c,e cos a c,e = < z £a,e " z ha,e>/ £ c ^ 

The tension force is resolved into helicopter body axes by the expression 
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The moments acting on the helicopter due to the cable tension are computed 
according to 


t,h Z t,h^a,h ~ Y t,h z a,h 
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The tension force is resolved into load body axes by the expression 
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and the moments acting on the load due to the cable tension are given by 


L t,£ = z t,£ y a,£ “ Y t,£ z a,£ 

M t f £ = x t,£ z a,£ " Z t,£ x a,£ / 
N t,£ = Y t,£ x a,£ " X t,£ y a,£ j 


(77) 


Equations of Motion 

The equations of motion for both the helicopter and the external sling 
load are developed in body axes with respect to a flat, nonrotating Earth. It 

is assumed for convenience that each body is rigid and that the x^-z^ plane 

and the x^-z^ plane are planes of mass symmetry and that gyroscopic effects 
of engines are negligible. The equations of motion for the helicopter are 
developed first. 

The translational motion equations for the helicopter are given in terms 
of body-axes components of angular velocity, translational velocity, accelera- 
tions, and components of gravity as 

”cg,h = r h v cg,h ' ^h w cg,h + a x,h ~ 9 sin 9 h 

^cg,h = Ph w cg,h - r h u cg,h + a y,h + 9 cos 6 h sin cf> h \ (78 
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The summations include all external forces: main- and tail-rotor forces , fuse- 

lage aerodynamic forces, and cable suspension forces. The last terms on the 
right-hand side of equations (78) give the acceleration components due to 
gravity. The inertial velocity of the helicopter is given in body axes by 
integrating equations (78) as 
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h dt + W cg,h (0) 
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The equations of rotational motion are used in the following form: 

Ph = [ ZL h - (I zz,h - I yy ,h ) ^h r h + I xz,hPh% 

+ f N h " ^ I yy ,h “ I xx,h^h^h ” I xz,h c ^h r 3q]^ [ /( I xx,h “ I xz,h/ I zz,h) 

zz,hj / 

qh = |^ M h ” ^xx,h “ *zz,h^h r h + I xz,h( r h ” ^h )J / X yy ,h 


> (81) 


*h [^ N h ■ I yy,h I xx,h^h < ^h + I xz,h^h c ^h r h^^ I zz,h 


J 


The body-axes components of angular velocity are determined by integrating 
equations (81) as 



\ p h dt + p h (0) 

q h = ’ 

^ q h dt + q h (0) 
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(82) 


The helicopter Euler angles shown in figure 3(a) are determined by 
integrating the following differential equations: 
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= ( % sin + r h cos < f ) h ) / cos 9 h ^ 

9 h = % cos " r h sin } (83) 

= Ph + (g h sin + r h cos V / tan 0 h > 

The inertial velocity of the helicopter c.g. expressed in body axes is 
given in Earth coordinates by the transformation 
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The Earth-axes coordinates of the helicopter c.g, are determined by integrating 
equation (84) to obtain 

x, = \ u, dt + x, (0) 

hcg,e J ncg,e hcg,e 

^hcg,e = J v hcg,e dt + / (85) 

z hcg,e = J w hcg,e dt + z hcg,e (0) ^ 

The equations of motion for the external load may be obtained by changing 
all the subscripts h to Z in equations (78) to (85) . The external- load 
equations of motion thus obtained may then be solved together with the equa- 
tions of motion for the helicopter to obtain load motions. 


APPLICATION TO U.S. ARMY CH-54 HELICOPTER AND CARGO CONTAINER 

Rotor and fuselage data for the U.S. Army CH-54 helicopter were obtained 
from unpublished Sikorsky Aircraft data, flight tests at Langley Research 
Center, and reference 9. The fuselage and main- and tail-rotor data used in 
the simulation are listed in table I. The fuselage wind-tunnel data are given 
in figures 5 to 11. These curves are entered with local fuselage angle of 
attack o if£, wind-tunnel yaw (sideslip) angle local incidence at the 

tail i t , and dynamic pressure q^ as determined from equations (53) , (56) , 
(54), and (50), respectively. Values of e kt and e^ f of 1.8 and 0.5, 
respectively, were used in these equations to determine fuselage lift, side 
force, rolling moment, pitching moment, and yawing moment. 
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The fuselage lift is calculated by using the data from figures 5 and 6 
and the following equation: 


L h “ 


AL ih A£ 2 h\-, 
— — + ~ — hh 

, q h %/ 


( 86 ) 


The fuselage wind-tunnel pitching moment is calculated by using the data 
from figures 9 and 10 and the following equation: 
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The fuselage drag is calculated according to the formula obtained from 
Sikorsky Aircraft 


D h = ( 7.25 + 2.4a fZ + 42.9a 2 fZ + 45.6^ t )i h 


( 88 ) 


where a f ^ and are in radians. 

The fuselage moments due to body angular rates are computed according to 
the formulas obtained from unpublished Sikorsky Aircraft data as 


L d , h 9 5 * 6r h V a s , h 

M d,h “ ' 218q h V as f h ^ 
Nd,h = -322r h V aS/h ^ 


(89) 


where V ag ^ is computed by using equation (51) , and r^ and q^ are in 
radians per second. Note that the rolling moment h is a function of r^ 
not p^. The fuselage aerodynamic forces and moments thus determined are used 
in equations (58). It is noted that in the transformation of aerodynamic 
forces from wind-tunnel to body axes given by equation (57) , the sideslip 
angle 3f was inadvertently set to zero. The change in handling qualities 
caused by this error is considered to be negligible at cruise conditions and 
nonexistent at hover. 

Mass data, cable attachment point coordinates, ground contact coordinates 
and parameters, and nominal cable length used for a 2.4-m by 2.4-m by 6.1-m 
cargo container are given in table II. Aerodynamic data used for the cargo 
container are derived below from data given in references 10 and 11. Since 
large angles of sideslip and angles of attack were expected during cruise and 
hover flight, the wind-tunnel data - which ranged from -5° to 45° in pitch and 
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from 0° to 95° in yaw - had to be modified to allow ±180° angle-of-attack varia- 
tion and ±90° angle-of-sideslip variation* This was accomplished by fitting 
trigonometric functions to the existing data, as shown in equations (92) 
and (94). 

The load components of free-stream velocity u_ 0 , v„ Q , and w„ Q are 

I X/ aLS r Xf do / X/ 

obtained by adding steady-wind, gust, and inertial components calculated by 
using the load version of equations (4), (10), and (80) as follows: 


U as,£ ~ u cg,£ + u gust,£ + U wind,£ 
v as,£ = v cg,£ + v gust , £ + v wind,£ / 
w as,£ ~ w cg,£ + w gust,£ + w wind,£ 


(90) 


Then the load angles of attack and sideslip and the dynamic pressure are calcu- 
lated from the following expressions: 
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It was found in simulation that during hover with no winds, divergent oscilla- 
tions of load angular motions would occur. This was due to the fact that in 
the zero-airspeed condition, load oscillations due to cable stretch were essen- 
tially undamped and the slight phase shift due to the numerical integration of 
the equations of motion caused a numerical instability. This instability was 
eliminated by adding a constant "rotor downwash" value to the z^-component of 
velocity as follows: 


w n = w 0 
as,£ as,£ 


9.14 


The cargo container lift, drag, and side forces in wind axes at the load geo- 
metric center are calculated from expressions derived from the data in refer- 
ence 10 as 
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L£ = (6.5 sin 2a^ cos 

= [20.9 - 7.66(1 + cos 2a^ cos B£)J q£ 
= (-7.9 sin 2g £ cos 2a^)q£ 




( 92 ) 


These forces are expressed in load body axes by using the transformation 
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The cargo container aerodynamic moments are given below and are based on the 
assumption that the static pitching and rolling moments are negligible in com- 
parison with the pitching and rolling moments caused by the suspension system. 
The static yawing moment was derived from the data in reference 10 as 




; % ~ ^-7.73 sin 23£ + 4.47 sin 43^ cos IK 


(94) 


The aerodynamic moments due to rolling, pitching, and yawing angular 
velocities are assumed to be proportional to the product of airspeed and angu- 
lar velocity, as in the case of the fuselage moments given by equations (89) , 
and are given by 

L d,H = -89V as,£ p Jl 

M d,* = - 89V a S/ *5£ > (95 

N d,& = “ 89V as,£ r 2, > 


where p£, q£, and r £ are in radians per second, 

moments acting on the load are thus given by 


The total aerodynamic 
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SIMULATION DESCRIPTION 


Computer Hardware 

The mathematical model has been programed in FORTRAN IV for operation on 
the Langley real-time simulation system using the Control Data CYBER 175 
digital computer system. The program accepts inputs from the simulation cock- 
pit through ADC's, and outputs voltages to the simulator through DAC's. The 
program flow is controlled by an operator through use of a simulation control 
console (fig. 12) . 

The simulation control console consists of white indicator lights, red 
indicator lights, function sense switches, mode control switches, a data entry 
keyboard, digital decimal display unit, and potentiometers. The white indicator 
lights are used to indicate program status or diagnostics. The red indicator 
lights are used to indicate program diagnostics. The function sense switches 
are used to select program options. The mode control switches shown in fig- 
ures 13 and 14(a) are used to control the running of the RTS computer program. 
Each switch (mode) is briefly described as to its use (mode nominally active 
when switch depressed) : 

OPER (OPERATE) - begins updating time and integrating the differential 
equations 

HOLD - holds integrated variables at previous value 

RESET - initializes program at Time = 0 

IDLE - idles the computer (no computations) 

CHANGE - changes program variable to the new value entered on the data 
entry keyboard 


SCAN - scans through tables and displays values on the digital decimal 
display unit 

RELEASE - releases CHANGE and SCAN modes 
ERASE - erases real-time disk file 

TERM (TERMINATE) - terminates program at simulation control console and 
transfers control to the Tektronix terminal 

READ - loads read overlay 

PRINT - loads print overlay 

RELEASE - releases ERASE, TERM, READ, and PRINT modes 

The data entry keyboard shown in figures 13 and 14(b) is used to input new 
values for program variables. The keyboard is used in conjunction with the 
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digital decimal display shown in figures 13 and 14(c), Any program variable 
can be changed with these two units by following a simple procedure. 

The potentiometers are used to input variables (through ADC’s), They are 
mainly used for checkout. 

The Tektronix terminal is used to communicate interactively with the 
CYBER 175, When the RTS program is running, the terminal is used as an output 
device, for example, for error messages. 


Cockpit 


The cockpit, which is linked to the CYBER 175 computer for the simulation, 
is shown in figures 15 to 18. It has simulated instrumentation, pilot controls, 
and a visual landing display system to simulate a helicopter cockpit. 

The cockpit instrument panel is shown in figure 15. Figure 16 shows the 
controls by which pilot inputs are fed into the CYBER 175 computer. They 
include the standard helicopter controls: cyclic stick with trim release 

button and "coolie hat" trim control, collective stick, and antitorque pedals. 

The cyclic stick is operated by a three-axis hydraulic control loader 
which is controlled by the analog computer shown in figure 19. Stick travel at 
the center of the hand grip is limited by mechanical stops to ±12.7 cm in the 
lateral direction and ±14.0 cm in the longitudinal direction. The stick 
dynamics are modeled on the analog computer as second-order systems with the 
damping chosen by the pilot. The stick force gradient used is 4.4 N/cm in the 
longitudinal axis and 8.9 N/cm in the lateral axis and was obtained from ref- 
erence 9. The force the pilot applies to the stick is opposed by the hydraulic 
control loader and measured by force transducers. The force transducer signals 
are fed into the analog computer and the stick acceleration is calculated. The 
acceleration signal is processed together with the actual stick position and a 
position error signal is formed. This error signal is then sent to the control 
loader servo which moves the stick. The new stick position is fed through an 
ADC to the CYBER 175 computer as the pilot’s cyclic stick input. The lateral 
position of the stick is the input x^ at in equations (12) and the longitudinal 
position is the input x lon in equations (12) . 

In the trimming mode, the CYBER 175 computer uses the trim algorithms dis- 
cussed subsequently in the section "Trim Calculations" to obtain the trim posi- 
tion of the cyclic stick. This position is fed through a DAC to the control 
loader analog computer as a command to drive the stick to the trim position. 

Activating the trim release button on the cyclic stick or stick trim 
toggle switch shown in figures 16 and 18 removes the trim force. The stick 
position at the instant the pilot releases the button or turns the stick trim 
toggle on becomes the new trim position. The trim position may also be moved 
by the four -position coolie hat mentioned previously. By moving the coolie hat 
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in a certain direction , the stick Vs zero force position is moved in that direc- 
tion* The trim release button, the stick trim toggle, and the coolie hat are 
discrete inputs to the CYBER 175 computer* 

The analog output from the collective stick, the quantity x cq1 in equa- 
tions (12), is fed through an ADC to the CYBER 175 computer. The stick travel 
is limited by mechanical stops to 35.6 cm. A friction collar can be adjusted 
to vary the stick friction. The trim algorithm calculates the collective-stick 
trim position. The pilot must then place the collective stick in its trim 
position, which is evidenced by the turning off of the collective trim light 
shown in figure 18. 

The pedals are driven by the same three-axis hydraulic control loader sys- 
tem which drives the cyclic stick; however, the response characteristics are 
modeled on the CYBER 175 computer because the frequency response of the pedals 
is not as high as that of the cyclic stick. The pedal travel is limited by 
mechanical stops to approximately 10 cm. The pedal dynamics are second order 
with a force gradient of 62.2 N/cm, as obtained from reference 9. The pedal 
damping was chosen by the pilot. The operation of the pedal control loader is 
identical with that of the cyclic stick. Microswitches located on each pedal 
have the same function as the trim release push button located on the cyclic 
stick. In the trimming mode, the trim algorithm calculates the trim position 
for the pedals. These position commands are fed through DAC's to drive the 
pedals to their trim position. 


The remaining pilot controls in the cockpit are discrete inputs to the 
CYBER 175 computer. They are 

Winch control - three-position toggle switch to lengthen or shorten the 
load cable 

Mode control buttons - shown in figure 18; select one of the following 
modes of operation: 

RESET resets to initial flight conditions 

HOLD holds present flight conditions 

OPERATE starts simulation 

AFCS selector buttons - shown in figure 18; select the following options: 


AFCS 

automatic flight control system on 

YAW 

heading hold on 

ALT 

altitude hold on 

HOVER 

hover hold on 
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Stick trim switch - two-position toggle switch to turn cyclic-stick-force 
trim system on and off; shown in figure 18 

Load release button - located on the cyclic stick; when pressed, external 
load is released 

Several additional features are included in the cockpit. 

The pilot views a visual terrain scene provided by the visual landing dis- 
play system (VLDS) , and a computer-generated load/landing zone scene. Pilots 1 
comments during checkout indicated that a separate and independent computer- 
generated video display of the sling load and landing zone as viewed by a down- 
looking TV camera attached to the helicopter at the pilot's location is unaccept- 
able. Many hours of simulation experience with NASA pilots and pilots from 
industry have shown that an acceptable simulation sling load display can be 
made by electronically combining the visual terrain scene with the computer- 
generated scene. The resulting image is then displayed on the video monitor 
placed in the normal eye-level position. Figure 17 shows the virtual image 
lens system through which the pilot views the color monitor which displays the 
visual terrain scene electronically mixed with the down-looking load scene. 

The virtual image lens causes the pilot's eyes to be focused at infinity. Also, 
shown in figure 17 is a small black and white monitor which shows the same dis- 
plays as the pilot's color monitor. The black and white monitor is used by the 
researcher sitting in the right seat of the cockpit. 

An audio generator creates cockpit noise, with one component having a 
frequency equal to 6 times rotor speed {in rpm) and the other component having 
white noise with magnitude proportional to the quantity (V h + 24.4) m/sec 
to represent wind noise. These audio signals are output on a speaker in the 
cockpit with volume level based on pilot comments. 


Simulation So f tware 

The simulation software is written in FORTRAN, with some assembly language. 
It is broken down into primary and secondary overlays with many subroutines and 
general-purpose functions. 

Figure 20 shows the central memory layout with the various overlays. The 
large block labeled (0,0) represents the base overlay, which requires 40 000 3 
storage locations. It contains all variables which must be communicated 
between overlays, subroutines required by more than one overlay, and some real- 
time system software. 

The primary level overlays (1,0) , (2,0), (3,0), and (4,0) share the same 
memory locations. Their functions are necessarily mutually exclusive. Over- 
lay ( 1 , 0 ) initializes program variables and initializes the real-time system. 
Overlay (2,0) prints trim sheets or time histories. Overlay (3,0) reads in 
data for various loads. Overlay (4,0) is the main overlay and contains the 
real-time loop. It is further described below and in figure 21(a). Secondary 
level overlays (4,1) to (4,10) are optionally executed by overlay (4,0). 
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The secondary overlays share the same memory locations. Overlay (4,1) 
executes static checks; (4,2) contains the load/landing zone display equations; 
(4,3) contains the trim algorithm and logic; (4,4) contains preprocessing equa- 
tions, i.e. , calculations which can be done out of the real-time loop for a 
given run; (4,5) checks and prints out function data; (4,6) calculates sta- 
bility derivatives; (4,7) prints error messages to the console operator; 

(4,10) calculates linearization derivatives. Of these eight functions , only 
the load/landing zone display equations are needed in the real-time loop. The 
others are brought in as necessary between runs. 

The real-time overlay (overlay (4,0) ) is the main loop, with all the other 
overlays supporting it. A general flow diagram of the main loop is shown in 
figure 21(a). This real-time loop consists of three separate loops: reset, 
hold, and operate. The reset loop is cycled through at Time = 0. The hold 
loop is selected at any time to hold all simulation variables fixed. The 
operate loop begins the simulation, calculating and integrating the equations 
of motion and updating the independent variable time in synchronization with 
real time, with a step size At of 1/32 sec. This rate was selected as the 
largest time interval allowed by the real-time system which gives accurate 
results and updates cockpit instruments with no visible jumps. 

The block in figure 21(a) labeled (A) represents the section of the reset 
loop in which various secondary overlays can be selected and executed. For 
example, when the console operator depresses the trim button, overlay (4,4) is 
executed to calculate preprocessing equations, and then overlay (4,3) is 
loaded to execute the trim algorithm. Overlay (4,3) will stay in memory 
until another option is selected. 

Blocks (B) to (E) are cycled through in both hold and reset modes. 

Block (B) represents the sampling of 5 analog signals (via ADC's) and 27 dis- 
crete inputs from the cockpit. These voltages are scaled and the resulting 
parameters (e.g., stick inputs) saved for later use. 

Block (C) represents all the calculations necessary to compute the equa- 
tions of motion for the helicopter and, optionally, for the load. This block 
is broken down in figure 21(b), where each block refers to a portion of the 
previously developed mathematical model. Some of the general-purpose algo- 
rithms used in this section are a convolution integration scheme used to repre- 
sent first- and second-order filters, a linear interpolation scheme for func- 
tion data table look-ups for functions of one or two variables, and an Euler 
integrator for use when speed is desirable and accuracy not critical. 

Block (D) represents the calculations for the VLDS and Adage displays. 

The VLDS drive equations are computed, 19 analog signals are output, and 
11 analog error signals are fed back. The Adage display equations are com- 
puted, and the data are transmitted in digital form to the Adage Graphics 
Terminal (AGT 130) . Further descriptions of the VLDS and Adage displays are 
given below. 
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Block (E) represents the scaling of program variables and the outputting 
of those variables via DAC's and discrete channels. Twenty-nine voltages and 
12 discrete outputs are trunked to the simulator cockpit. 

Blocks (F) and (G) are only cycled through in the operate mode. Data for 
the time-history printout are recorded in block (F) , every specified number of 
iterations. These data can be printed later at the researcher's option. 

Block (G) represents the integration of the equations of motion. The 
real-time loop is designed to be as fast as possible, while maintaining good 
accuracy. The Adams-Bashforth second-order (AB-2) , one-pass scheme is very 
fast, and is accurate enough for this particular simulation, as was shown in 
the comparison of time-history plots with the independent check program. The 
AB-2 integration scheme is fully described in reference 12. 

The complete real-time loop requires a minimum of 2.5 msec (in reset) and 
a maximum of 6.2 msec (in operate). The 2.5-msec figure includes only the heli- 
copter with AFCS , whereas the 6.2-msec figure includes the 2.4-m by 2.4-m by 
6.1-m load, winds, gust, and so forth. Table III shows a complete breakdown of 
the required time . 

Visual Landing Display System 

The visual landing display system (VLDS) consists of a fixed, colored- 
terrain board system and a movable camera transport and is designed for use 
with a monitor and virtual image lens system for displaying an 11 out- the -window” 
scene in a simulated cockpit. A brief description follows. A more detailed 
description is given in reference 13. 

The 7.3-m by 18.3-m terrain model board of the VLDS includes two airports 
and surrounding terrain, one at 750/1 scale and the other at 1500/1 scale, and 
is shown in figure 22. There are a total of five paved runways, from 0.6 km 
to 3.5 km in length. A helipad is located on the 750/1 airport and is shown 
in figure 23. It consists of a Maltese cross with a 45-m by 45-m border. The 
terrain is generally flat, and provision is made for variable visibility, 
variable cloud-base heights, and day, dusk, and night scenes. 

The camera system has a field of view 48° wide and 36° high and uses a 
525-line color video raster system. This system provides nominal resolution 
on the order of 9 minutes of arc. The pilot's eye position and the orientation 
of the pilot's line of sight out the forward window are both calculated with 
respect to the simulated runway. These positions and orientation angles (and 
their rates) are used to drive the camera system. The dynamic characteristics 
of the VLDS transport system are given in table IV. 


Load/Landing Zone Visual Display 

The simulation is intended to be used by pilots to evaluate sling load 
stabilization systems, and since sling load piloting is a visual task, a 
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realistic display of load motions and landing zone position must be available 
in the simulated cockpit. Preliminary work with pilots defined the format of 
the simulation display as follows. The display should show in perspective the 
three-dimensional outline of the external load, the cable (s) , and the Earth- 
fixed features such as the landing zone and lead-in distance marks. The center 
line of the display should be located at the cable attachment point on the heli- 
copter and the view should be stabilized in pitch and roll; that is, the view 
should always be vertically oriented. Load altitude above ground level should 
be presented in analog form by a pointer and altitude scale, and helicopter 
pitch and roll attitude bars should be generated to provide the pilot precise 
attitude information needed while hovering. A sketch of the desired load/ 
landing zone display is given in figure 24. 

The load/landing zone computer graphic display is defined by vector end 
points. These end points are computed by using the formulas developed in equa- 
tions (97) to ( 10 1 ) . 

The coordinates of the load c.g. are given with respect to helicopter 
body axes by 


rx £cg,h"j 

r 

1 

/ y £cg,h 


v z £cg,h v 

V 


x £cg,e x hcg,e^ 


z £cg,e z hcg,e J 


(97) 


The coordinates of any point p on the load expressed in helicopter body axes 
are thus given by 




x £cg,h 


r v "\ 

x £p , a 

< Y £p,h 

1 

>-< 

y £cg,h 

J 

) + c h c u T < 

y £p, £ 

\Hp f hJ 


v z £cg,hJ 


<. z £p ,£> 


(98) 


The coordinates of the helicopter viewpoint expressed in Earth axes are given 
by 


r x 

x v,e 

y v,e 

)- < 

A 

hcg,e 

^hcg,e 

| 

> + c h T < 

( X 

v,n 

y v,h 

^- Z v ,e j 


< z hcg,ej 


^ z v,h > 


(99) 


45 



where x , , y , , z , are the coordinates of the viewpoint in helicopter 
v,n v,n v,n 

axes. Thus, the coordinates of any point p on the Earth with respect to the 
viewpoint are given in helicopter axes by 


fx ^ 


r x 

- X 

zp,h 


zp,e 

v,e 

y zp,h 


Y zp,e 

^v,e 

v. z zp,hj 


< z zp r e 

z v,e^ 


(100) 


where x„^ . y„„ z„_ . are the coordinates of any point p in Earth axes, 

zp , e zp , e zp , c 

Since the load/landing zone display will be shown on a screen, these points 
must be expressed in the coordinate system of the screen. The following equa- 
tions are used to transform these points to screen coordinates and to add 
perspective : 


x screen x p,h \ p,h 

\ ( 101 ) 

^screen “ ^p,h / z p,h ^ 

where x^ y^ ^ are the coordinates at any point p in helicopter 

axes. The helicopter attitude bars are generated directly in screen coordi- 
nates. The data defined in screen coordinates must be clipped at the screen 
boundary. The clipping algorithm employed is described in reference 14. 

The load/landing zone display requirements are met by the Adage Graphics 
Terminal in conjunction with the CYBER 175 computer. The Adage is an inde- 
pendent, digital, graphics computer with an operator console shown in figure 25. 
Computations required for the display are performed on the CYBER 175 computer 
and data are transmitted to the Adage, which generates the display. Sample 
displays are shown in figures 26 and 27. 


Trim Calculations 

In order for the simulation to start in an unaccelerated flight condition, 
it is required that an algorithm that renders certain mathematical model time 
derivatives zero be available for determining the helicopter attitude and con- 
trol positions and the load attitude and position with respect to the heli- 
copter. The algorithm used is based on the so-called method of secants and is 
described in appendix D of reference 15. The algorithm determines the values 
of the following independent variables: x lon , x i at ' x col' x ped f a f' ^h' ®h' 

v m' v t' 9 0t' a i' H’ Q V x £cg,e' y £cg,e' Hcg,e so that the following dependent 
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variables are approximately zero: u h , v cg ^ h , w h , p h , q h , 

v m' ^t' ®0t' u cg,£' v cg,Z f w cg,Z r $2.' r Jt* 


h' w hcg , e ' 


For the simulation of the CH-54 helicopter and cargo container the algo- 
rithm trims at any specified airspeed from -20 to 100 knots in approximately 
1 sec. Values of the major helicopter variables in trim at 0.1, 30, 60, and 
90 knots are given in table V. The use of 0.1 knot was to prevent singularities 
that would occur at exactly zero airspeed. 


VERIFICATION AND VALIDATION 
Simulation Verification 

The simulation software verification was accomplished with an independent 
check program. The check program was written in FORTRAN and was designed to 
run in a batch environment. It was developed from the mathematical model inde- 
pendently of the real-time program. As the independent check program is not 
required to run in real time, a more accurate integration algorithm than AB-2 
was selected for this use. This procedure provided a verification of both the 
programing of the mathematical model and the accuracy of the numerical solu- 
tion. A Runge-Kutta seventh-order, 13-pass integration algorithm with vari- 
able step size was used in the independent check program. This algorithm and 
its accuracy are described in reference 16. 

The check case selected consisted of a 30-sec run, during which the four 
control inputs were disturbed independently with a half-cycle sine-wave forcing 
function. The AFCS was off for the first 3 sec, and was cut on and off in 
3-sec intervals during the run. This allowed for large, but reasonable, per- 
turbations in the state variables. The check case was run for the helicopter 
alone at hover and at 60 knots and for the helicopter with sling load at hover 
and at 60 knots. The time histories of the state variables showed excellent 
agreement with those from the real-time program, and further analysis proved 
that the negligible differences were due to integration errors. These errors 
were small enough to verify that the AB-2 integration scheme was sufficiently 
accurate for the simulation. 

Hardware verification consists of a static check performed before each 
real-time session. The instruments are checked visually for conformity with 
static check values. All the control inputs are deflected to their maximum 
and minimum positions while the console operator monitors the voltages. 

Finally, all the discrete inputs are activated while the console operator 
monitors them. 


Mathematical Model Validation 

The helicopter portion of the previously described mathematical model was 
validated by comparison of simulated time histories with flight data, analysis 
of eigenvalues and eigenvectors, and by pilot evaluation. 
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Eleven simulations were made for comparison with flight time histories, 
and the results from two of these are shown in figures 28 and 29. In all 
cases the simulated helicopter mass characteristics were set equal to those of 
the helicopter that generated the flight time histories. (These flights were 
conducted at Langley in 1972 in an unrelated sling load project reported in 
ref. 17.) Table I gives the mass characteristics and other data used in the 
simulation. The AFCS was off in all cases. When reviewing figures 28 and 29, 
the reader should keep in mind that the basic helicopter is unstable about each 
axis, so that any perturbation from trim will result in divergent roll, pitch, 
and yaw motion because of the dynamic cross coupling. This divergence tends to 
amplify differences between flight and simulation as time increases. It was 
the manufacturer's opinion that about 2 or 3 sec of close agreement between 
simulation and flight data would be all that should be expected. 

Figure 28 shows the data for a run in which the helicopter alone was 
hovering and the pilot made a longitudinal cyclic step input of about 3 cm for 
approximately 2 sec while holding the other controls approximately fixed. The 
simulation was trimmed at hover and the control position changes made in flight 
were duplicated in the simulation by using the recorded control input data. The 
upper two plots in figure 28 show the longitudinal and lateral cyclic stick 
motion for the first 5 sec of this case. The pedals moved less than 0.05 cm 
during this run and the collective stick moved less than 0.09 cm at the grip. 

The differences in Euler angles and body rates between flight and simulation at 
time zero noted in the figure were due to trim and flight data bias and should 
be disregarded throughout the run when comparing these data. Except for the 
slight roll acceleration reversal noted in the figure at about 3 sec, the data 
for attitude and angular rates show very close agreement. The helicopter alti- 
tude data agreed closely between simulation and flight and showed a decrease of 
about 1.5 m during the 5-sec period. 

Figure 29 shows the data for a run in which the helicopter was flying 
straight and level at approximately 38 knots and the pilot made a longitudinal 
cyclic step input of about 2 cm for several seconds followed by a pedal input 
of about 1 cm and a lateral cyclic input of approximately 2 cm. The same trim 
airspeed was used in the simulation, and the control position changes made in 
flight were duplicated in the simulation as before. The upper three plots in 
figure 29 show the time histories of the flight and simulation control inputs 
for the first 7 sec of this run. The collective stick was not moved during the 
7-sec period shown. The time histories of helicopter altitude and airspeed are 
given in the fourth and fifth plots, where good agreement is noted between 
flight and simulation. As for the hover run, the differences in Euler angles 
and body angular rates at time zero should be disregarded throughout the run 
when comparing the data. Here the agreement in pitch and yaw data is considered 
excellent, whereas the agreement in roll data is only fair. The remaining nine 
comparison cases between flight data and simulation data showed excellent agree- 
ment for the first few seconds in each case. In some cases the agreement 
remained good for several more seconds; in others the agreement degraded 
slightly, as in the examples- shown in figures 28 and 29. In general, the 
agreement between simulation data and flight data was believed to indicate that 
the simulation mathematical model is adequate for studying helicopter sling load 
con tr o 1 systems. 
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The helicopter was trimmed at 0.1, 30, 60, and 90 knots and the lineariza- 
tion algorithm described in reference 19 was applied to determine the linear 
system derivative matrices A and B defined in reference 19 and the system 
eigenvalues and eigenvectors . In order that the time-varying rotor inflow 
ratios V m and \> t and the tail-rotor effective pitch angle 0Q t would not 
appear as variables in the linear system, 30 iterations with the V m , V t , 
and. 0Q t differential equations (eqs. (26) and (47) ) were performed after each 
state variable perturbation was made so that these three variables would reach 
semi -steady- state before the linear system derivatives were calculated. The A 
and B linear system matrices are given in table VI for the four airspeeds. 
Unpublished derivatives obtained from Sikorsky Aircraft, derivatives calculated 
by methods given in reference 7, and derivatives given in reference 18 were 
compared with those in the table and generally good agreement was found. 

The eigenvalues and eigenvectors for the system A matrix were calculated 
and are also given in table VI. The upper half-plane of the loci of eigen- 
values at the four airspeeds is shown in figure 30. The eigenvectors were used 
to identify the eigenvalues with respect to the helicopter modes of motion. 

The identification is noted on the figure. The real roots ranging from -0.85 
to -1.0 are associated with rolling velocity and vertical velocity. The zero 
at the origin corresponds to the heading angle The root loci were com- 

pared with unpublished Sikorsky Aircraft data for a CH-53 helicopter, which is 
similar to the CH-54, and with data at 60 knots given in reference 18. In gen- 
eral, the loci trends were found to be similar but the agreement in magnitudes 
was only fair. 

The comparison of the linear system matrices and the eigenvalues and 
eigenvectors obtained from the nonlinear mathematical model with data from other 
sources indicates that the mathematical model does represent the dynamics of a 
CH-54 helicopter fairly well. 

The sling load portion of the mathematical model was validated by compar- 
ing measured frequencies of simulated pendulum, rocking, and bounce motions 
with theoretical values and by pilot evaluation. The theoretical and measured 
frequencies of these modes are given in table VII. In each case the theoreti- 
cal values are based on the assumptions that the helicopter attachment point is 
fixed in space and that the cable is inelastic except in the bounce mode. In 
each mode the measured frequencies are slightly higher than the theoretical 
values. This discrepancy is believed to be due to the motion of the helicopter 
attachment point and the interaction between the bounce mode and the rocking 
modes. The agreement is considered adequate. 

As mentioned earlier, the suspension system cable spring rate K was 
critical in terms of numerical stability. Trial and error showed that the 
spring rate had to be selected so that the frequencies associated with the 
cable stretch would be no higher than 2 Hz. A value of K gc = 1.8 x 10 5 N/m 
was used with the 45 36 -kg load, giving a theoretical bounce natural frequency 
of 1 Hz. 
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Pilot's Comments 


The project test pilot acquired experience in the Army CH-54 helicopter 
during the flight tests reported in reference 17. This pilot also received 
additional flight training with long-line sling loads in connection with the 
present report. Numerous simulated flights both with and without a simulated 
sling load attached were performed for the purpose of validating the simulation 
mathematical model. 

Several combinations of simulated cyclic stick and pedal force gradients 
and dynamic characteristics were considered. The set used in this report repre- 
sent a usable set for the simulation. The VLDS and load/landing zone visual 
display system described earlier is considered to be marginally adequate for 
the simulation of helicopter sling load operations. The lack of peripheral 
visual cues in the horizontal and vertical planes makes hover and transitions 
to hover especially difficult. The computer-generated load/landing zone dis- 
play and the VLDS also lack the depth and natural cues that are so useful in 
actual flight. 

The simulated helicopter alone with and without winds and turbulence seemed 
like a large helicopter in all flight regimes. It seemed underdamped irr yaw 
motions when approaching hover, and its vertical motions seemed much more sen- 
sitive to collective movements than those for the actual helicopter. The simu- 
lated helicopter seemed slower than the actual helicopter to respond to cyclic 
step and pulse inputs. Autorotation was simulated and seemed realistic to 
touchdown. 

The simulation with the external load seemed realistic in spite of the 
visual display shortcomings. Simulated flights with the external load were 
performed over the entire speed range and the ensuing load motions were observed 
on the load/landing zone display. In all cases the load motions appeared real- 
istic. The helicopter motions definitely were affected by load motions, and 
hovering over a spot was very difficult, as it is in actual flight. All normal 
sling load operations were tried with and without winds and turbulence and were 
found to be possible to execute but with slightly higher work load than in 
actual flight. Simulator motion validation was not necessary since it was 
thought that for studies of sling load stabilization and control system com- 
parison, simulator motion would not be required. 

It is thought that the simulation described in the report can be used to 
compare various systems for stabilizing helicopter sling loads and improving 
helicopter sling load handling qualities. 


CONCLUDING REMARKS 

A generalized, full-flight envelope, real-time, piloted visual simulation 
of a single-rotor helicopter, suspension system, and external load is described 
and validated for the full flight envelope of the U.S. Army CH-54 helicopter 
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and cargo container as an example. The mathematical model described uses modi- 
fied nonlinear classical rotor theory for both the main rotor and tail rotor, 
nonlinear fuselage aerodynamics, an elastic suspension system, nonlinear load 
aerodynamics, and a load-ground contact model. The implementation of the 
mathematical model on a large digital computing system is described, and vali- 
dation of the simulation is discussed. The mathematical model is validated by 
(1) comparison of flight data with simulated data; (2) comparison of linearized 
system coefficient matrices, eigenvalues, and eigenvectors with calculated 
values, manufacturer's data, and data obtained from flight tests; and (3) by 
pilot evaluation. A visual landing display system that generates the pilot's 
forward-looking real-world display is discussed, and a special head- up, down- 
looking load/landing zone display is described. It was the test pilot's 
opinion that the simulation described in this report can be used to compare 
various systems for stabilizing helicopter sling loads and improving helicopter 
sling load handling qualities. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 23, 1978 


51 



REFERENCES 


1. Daniels, Glenn E . , ed. : Terrestrial Environment (Climatic) Criteria Guide- 

lines for Use in Aerospace Vehicle Development, 1973 Revision, NASA 
TM X-64757, 1973, 

2. Wilcock, T.; and Thorpe, Ann C. : Flight Simulation of a Wessex Helicopter - 

A Validation Exercise. C.P. No. 1299, British A.R.C. , 1974. 

3. Seckel, Edward; and Curtiss, H. C. , Jr. : Aerodynamic Characteristics of 

Helicopter Rotors. Rep. No. 659, Dep. Aerosp. Mech. Sci. , Princeton 
Univ. , Dec. 1963. 

4. Bailey, F. J. , Jr.: A Simplified Theoretical Method of Determining the 

Characteristics of a Lifting Rotor in Forward Flight. NACA Rep. 716, 

1941. 

5. Amer, Kenneth B. ; and Gustafson, F. B. : Charts for Estimation of 

Longitudinal-Stability Derivatives for a Helicopter Rotor in Forward 
Flight. NACA TN 2309, 1951. 

6. Gessow, Alfred; and Myers, Garry C. , Jr. : Aerodynamics of the Helicopter. 

Macmillan Co. , c.1952. (Republished 1967 by Frederick Ungar Pub. Co.) 

7. Seckel, Edward: Stability and Control of Airplanes and Helicopters. 

Academic Press, Inc., c.1964. 

8. Amer, Kenneth B. : Theory of Helicopter Damping in Pitch or Roll and a 

Comparison With Flight Measurements. NACA TN 2136, 1950. 

9. Bailes, Edward E. ; Diekmann, Vernon L. ; Watts, Joseph C. ; and Henderson, 

John C. : Instrument-Flight-Rules Capability Evaluation - CH-54B (TARHE) 

Helicopter. USAASTA Pro j . No. 71-01, U.S. Army, Dec. 1972. (Available 
from DDC as AD 908 656.) 

10. Laub, Georgene H.; and Kodani, Hifu M. : Wind Tunnel Investigation of 

Aerodynamic Characteristics of Scale Models of Three Rectangular Shaped 
Cargo Containers. NASA TM X-62,169, 1972. 

11. Liu, David T. : In-Flight Stabilization of Externally Slung Helicopter 

Loads. USAAMRDL Tech. Rep. 73-5, U.S. Army, 1973. 

12 . Wilson, John W. ; and Steinmetz, George C. (With Appendix A by Roland L. 
Bowles) : Analysis of Numerical Integration Techniques for Real-Time 

Digital Flight Simulation. NASA TN D-4900, 1968. 

13. Parrish, R. V.; Rollins, J. D. ; and Martin, Dennis J. , Jr.: Visual/Motion 

Simulation of CTOL Flare and Touchdown Comparing Data Obtained From Two 
Model Board Display Systems. AIAA Paper 76-1709, Apr. 1976. 


52 



14. Newman, William M. ; and Sproull, Robert F. : Principles of Interactive 

Computer Graphics, McGraw-Hill Book Co,, 1973, pp a 123-124, 

15. Houck, Jacob A.; Gibson, Lucille H. ; and Steinmetz, George G. : A Real-Time 

Digital Computer Program for the Simulation of a Single-Rotor Helicopter, 
NASA TM X-2872 , 1974, 

16. Fehlberg, Erwin: Classical Fifth-, Sixth-, Seventh-, and Eighth-Order 

Runge-Kutta Formulas With Step size Control. NASA TR R-287, 1968. 

17. DiCarlo, Daniel J. ; Kelley, Henry L. ; and Yenni , Kenneth R. : An Exploratory 

Flight Investigation of Helicopter Sling-Load Placements Using a Closed- 
Circuit Television as a Pilot Aid. NASA TN D-7776, 1974. 

18. Tomaine , Robert L. : Flight Data Identification of Six Degree-of-Freedom 

Stability and Control Derivatives of a Large "Crane" Type Helicopter. 

NASA TM X- 73 95 8, 1976. 

19. Dieudonne , James E. : Description of a Computer Program and Numerical Tech- 

nique for Developing Linear Perturbation Models From Nonlinear Systems 
Simulation. NASA TM-78710, 1978. 


53 



TABLE I.- VALUES OF PARAMETERS FOR CH-54 HELICOPTER 
[subscripts m and t denote main and tail rotor] 


x mr,h' 

m 

y mr,h' 

m 

z mr ,h' 

m 

x tr , h ' 

m 

Y tr,h' 

m 

z tr,h' 

m 

x wt,h' 

m 

y wt , h ' 

m 

z wt,h' 

m 

x ps,h' 

m 

-^ps,h' 

in 

Z ps,h' 

m 


FSCG, in 
WLCG, m 
BLCG, m 


x a,h' m * * 

y a ,h' m • * 

z a,h' m * * 
kg_m2 

*yy,h’ k 9- m2 
1 zz,h' kg - m2 
I xz,h' kg " m2 
e na/ m • • • 
t kg-m 

*W k< 3~ m 

® • * • 

T Am' sec ‘ 

6 3m' rad * 
0 lm' rad . 


. -0.33 

. 0 

-2.26 
. -13.74 

-0.84 
, - 2.22 
-0.51 
. 0 
-0.37 
5.77 
-0.53 
0.11 

. 8.20 

, 4.28 

0 

0.33 

0 

0.24 
. 39 800 

2.04 x 10 5 
1.78 x 10 5 
. 11 400 

. 0.610 
4.75 x 10 3 
715.0 
10.97 
0.20 
0 

. -0.183 

0.115 
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TABLE I.- Continued 


0 , rad * . « . „ 
sm' 

<j> , rad * . . » . 
T snr 

f rpm • • « • • • 

G gov' N~V (rad/sec) 
K dgov , N-m/ (rad/sec) 
K m/ N-m/ (rad/sec) • 
K gov' N ~ m/ (rad/sec) 
I pt r kg-m 2 .... 

W k ^ m2 - * * * 

T eng' sec 

GgQ, rad/rad . . . 

G Af' rad / rad • • • 

Won' rad / m ‘ * * 

G Axlat' rad/m • ‘ * 
rad/m .... 


G Q _, rad/ (rad/sec) 
tsq 

G Ap' rad /( rad / sec ) 
G 0tr' rad / (rad/sec) 

Gohm rad/rad . . . 


K c0' rad * 
K cl , rad/m 

K c2 , rad/m 

K c3 , rad/m 

K c4 , rad/m 


K c6 , rad/m 
K 7, rad/m 
K fe , N-m/N 

K sc' N/m ' 
m h , kg . . 

a t , per rad 


. . -0.0524 
. . 0 
. . 184.5 

. . 85 160 

1.32 x 10 5 
1.572 x io 6 
. * 833.3 

. . 4325 

. . 31 310 

. . 0.50 

. . 0.281 
. . -0.133 

. . 0.363 

. . 0.475 

. -0.00037 

. . 0.727 

. . -0.096 

. . 0.335 

. . 0.133 

. . 0.128 
. . 0.955 

, . 1.361 

. . -0.096 

. . 0.824 

. . 0.0494 

. . 3.64 

. . 1.09 

. . 0.0243 

1.8 x 10 5 
. . 13 610 

. . 5.73 
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TABLE I.- Concluded 


c t , m 


e t' 

m . . 

w 

kg-nr 

rt 

kg-m 

T At' 

sec 


T <53t' sec 
0 lt' rad 
R^, m . . 

CJ-j- . e * 

6 st' rad 
* s t' rad 
^2t r ra< ^ 


a m , per rad 
0 ) af rad/ sec 

p a ■■•••• 
e kt * * * • 
e kf 

a to • * • * 


0.92 

4 

0.343 

0.127 

13.88 

11.84 

0.20 

0.20 

-0.140 

2.44 

0.179 

0 

1.57 

0.78 

0.97 

6 

0.661 

5.73 

14.0 

1.0 

1.8 

0.5 

0 
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TABLE II.- VALUES OF PARAMETERS FOR 2.4-m BY 2.4-m BY 6.1-m 

CARGO CONTAINER 

I xx £/ kg-m 2 ........... 1124 

I yy l> kg-m 2 14 610 

I 2Z £, kg-m 2 14 610 

I xz£' k ^" m2 ‘ 0 

kg . 4536 

x a, £' m 0 

m 0 

z a, £' m -6.1 

x cl,£' m 3.05 

y c l , l ' m . . . , -1.22 

z cl / £' m .............. 1.22 

x c2,£' m 3.05 

y c2 ,V m 1.22 

Z c2 «•••••••••»*.••«.••••*••»»•».•» 1.22 

x c3,£' m * -3.05 

* 03 ,*' m * * 1.22 

z c3,Z' m ........ 1.22 

x c4 r Z' m -3.05 

y c4,Z' m ............. -1.22 

z c4,r m . . . 1.22 

V Q , m/sec ................... 0.305 

U f 0.25 

, N/m 45 000 

, N/ (m/sec) 7100 



TABLE III.- REQUIRED COMPUTER TIMES 


Calculation 

Computer 

time, msec. 

for mode - 

Reset 

Hold 

Operate 

Helicopter only 

2.5 

1.7 

2.0 

Cockpit interface 

.3 

.2 

.3 

Winds and gusts 

.0 

.1 

.2 

Visual landing 

.7 

.7 

.8 

display system 




Load/ suspension 

.9 

.9 

1.0 

system 




Load landing zone 

1.5 

1.4 

1.9 

visual display 




Total simulation 

5.9 

5.0 

6.2 


TABLE IV.- DYNAMIC CHARACTERISTICS OF VLDS 



.7 










TABLE V.- TRIMMED VALUES OF MAJOR VARIABLES AT FOUR AIRSPEEDS 
[Subscripts m and t denote main and tail rotor] 
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TABLE V.- Continued 


\7;=i V*n o 




Value at airspeed, knots 

, of 

- 








0.1 

30 

60 

90 

M mr,h' x- 1 ” 

• 


• 


• 

-9 39 

-3420 

-7734 


-8927 


N mr , h ' N " m 

• 


• 

• 

• 

1.20 x 10 5 

9.44 x 10 4 

7.54 x 

10 4 

8.37 x 

10 4 

X tr,h' N * 

• 


* 

• 

m 

-0.62 

-99.6 

-91.4 


-138.4 


Y tr,h- N ' 

• 


O 

• 

t 

8699 

6868 

5414 


5912 


Z tr,h' N * 

• 


• 


• 

0.17 

25.3 

19.4 


29.5 


L tr ,h' N " m 

• 

• 

• 

« 

• 

1.93 x 10 4 

1.53 x 10 4 

1.21 x 

10 4 

1.33 x 

10 4 

M tr , h ’ N “ m 

• 

• 

• 

• 

* 

-2288 

-721 

-254 


59.2 


N tr ,h' N " m 

• 


• 

• 

• 

-1.20 x 10 5 

-9.4 x 10 4 

-7.4 x 

10 4 

-8.09 x 

10 4 

X f,h' N • 

m 

• 

• 


• 

-0.010 

-1161 

-4401 


-9904 


Y f,h' N • 

• 

• 

9 


• 

0 

0 

0 


0 


Z f,h' N ‘ 

m 

• 

9 


o 

-0.018 

1046 

376.5 


865.3 


L f,h' N ~ m 






0 

0 

0 


0 


M f,h' N_m 

• 

9 

9 


♦ 

3227 ! 

4142 

7988 


8868 


N f,h' N - m 

• 

9 

* 


• 

-0.003 

-340.1 

-1252 


-2824 


f m • • 

• 


* 

• 

• 

30.5 

30.5 

30.5 


30.5 


P, kg/m 3 . 
q, N/m 2 

* 

• 

9 

• 

9 

• 

9 

9 

1.23 

0.0016 

1.23 

145.6 

1.23 

582.3 


1.23 

1310 


fyn' r P m • 

• 

• 


• 

9 

184.5 

184 .5 

184.5 


184.5 


C Tm * * * 

9 

• 



9 

0.00640 

0.00646 

0.00642 


0.00643 


P m . . . . 

• 

• 

• 

• 

9 

2.42 x 10 -4 

0.0726 

0.145 


0.216 


v m • • • • 

• 


• 

• 

9 

0.0566 

0.0388 

0.0217 


0.0146 


\n • * • • 

* 


• 

• 

9 

-0.057 

-0.041 

-0.031 


-0.042 


2 a m' N " m * 

• 

• 

• 

• 

• 

1.19 x 10 5 

| 

9.43 x 10 4 

| 

7.51 x 

10 4 

8.33 x 

10 4 
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TABLE V* - Concluded 



Value at airspeed, knots, of - 

variaDie 

0.1 

30 

60 

90 

V N 

1.33 x 10 5 

1.34 x lo 5 

1.33 x 10 5 

1.33 x 10 5 

H m» ** 

89.9. 

2381 

4198 

6628 

J m- « 

5.04 

1135 

1683 

2775 

rpm ...... 

835.6 

835.6 

835.6 

835.6 

C Tt 

0.00838 

0.00662 

0.00522 

0.00570 



2.41 x 10 -4 

0.0724 

0.145 

0.217 

v t 

0.0647 

0.0400 

0.0179 

0.0131 



-0.065 

-0.040 

-0.018 

-0.013 

2at' N " m 

2291 

1291 

724 

65 3 

T t / « 

8699 

6868 

5414 

5912 

«t' N 

0.624 

100.2 

92.2 

140.4 

J t , N ....... 

0.16 

22.4 

15.1 

17.2 

p h , 1/sec 2 

-3.6 x 10 -7 

-6.5 x 10~ 9 

3.8 x 10~ 8 

-1.4 x 10 -7 

q h , 1/sec 2 

1.2 x lo -6 1 

* -1.9 x 10~ 7 

2.3 x 10~ 8 

5.4 x 10 -8 

r h , 1/sec 2 

5.9 x 10 -7 

-1.2 x lo -8 

8.4 x 10~ 8 

-2.8 x 10 -7 

u^, m/sec 

2.1 x lo" 9 

2.3 x 10 -8 

2.1 x -lo -9 

2.4 x 10“ 9 

v^, m/sec 2 . . • * . 

-8.5 x lo -9 

2.4 x 10“ 12 

8.5 x 10“ 8 

3.3 x 10 -7 

Wft, m/sec 2 . • . . . 

1.5 x lo -9 

6.9 x lo -8 

1.8 x 10“ 9 

2.1 x 10 -8 

w hcg , e ' m / sec • * • 

-3.4 x 10 -12 

-1.8 x 10 -11 

-5.2 x lo -11 

-4.2 x 10“ 10 

0Q t , deg/sec 

-2 x io~5 

-2.5 x io~ 5 

-1 x 10“5 

0 

• 

Vrv, 

-1.5 x io~7 

0 

0 

0 

V t 

-1 x l(T 7 

0 

0 

0 
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TABLE VI*- HELICOPTER LINEAR SYSTEM MATRICES, EIGENVALUES 


AND EIGENVECTORS 


(a) Values at 0.1 knot 


lOClTY ■ 

.1 KNOTS 7 

POINT 

FORMULA 30 ITERATIONS 

INDEPENDENT VARIABLE 

TOTAL 

INCREMENTS 

.. "cg,h__ 

.30486-01 

*ic 

.17456-02 

— lcg t h 

.30486-01 


.17456-02 

”cg,h 

.30406-01 

8°m 

•17456-C2 

— 

.17456-02 

. 0 ct 

.17456-02 

q h 

.17456-02 



*h 

.17456-02 



_*h 

.17456-02 



6 h. 

.17456-02 




*17456-02 




4 MATRIX 



U cg , h 

V cg,h 

W cg,h 

P h 

q h 

r h 

*h 

6 h 


ycg,h 

-.1 3066-01 

-.73276-02 

-.83046-02 

-.16476-01 

.25276+00 

.17016-02 

.1824E-09 

-.98046+01 

-.14876-19 

Y c 8* h 

. 71216-02 

-.2995E-01 

-.51386-02 

.15196+00 

-.75436+00 

.19666+00 

.97926+01 

-.10566-01 

.14876-18 

w cg,h 

-.92046-02 

-.59046-02 

-.33376+00 

-.17686-01 

-.37896-01 

-.10126-01 

.47756+00 

.21656+00 

-.7667E-17 

p h 

.13616-01 

-.34446-01 

-.36846-02 

— • 7563E+00 

— .1262E+01 

.78116-01 

-.53676-06 

-.50826-20 

-.50826-20 

k 

.52726-02 

.25546-02 

-.62076-02 

.12546+00 

-.21706+00 

-.59946-02 

.3961E-07 

0. 

0. 


.10616-02 

.15656-01 

.83146-02 

-.22096-01 

-.67326-01 

-.24586+00 

.11786-05 

-.14916-18 

-.14916-18 

th 

0. 

0. 

0. 

. 1000E+01 

.10776-02 

-.22096-01 

0. 

0. 

0. 

9h 

0. 

C. 

0. 

0. 

.99886+00 

.40716-01 

0. 

0. 

0. 

*h 

0. 

0. 

C . 

0. 

-.48726-01 

.99916+00 

0. 

0. 

0. 


B MATRIX 



B IC 

A IC 

9 0m 

0 ct 

u cg> h 

.97546+01 

.72386-13 

-.20196+01 

-.54096-03 

?cg,h 

.28236-03 

•9767E+01 

-.15086+01 

•4165E+01 

”cg,h 

.52826+00 

-.76676-17 

-.90986+02 

.13786-03 

Ph 

.19656-03 

. 2015 E + 02 

-.74716+00 

,19416+01 

% 

-.38576+01 

8009E-10 

-.1 554E+01 

-.95766-01 

l h 

-.42456-03 

.11916+01 

.78576+01 

-.42606+01 

th 

0* 

0. 

0. 

0. 

9h 

0. 

0. 

0. 

0. 

*h 

0. 

0. 

0. 

0. 


REAL PARTS OF 

EIGENVALUES 








ROLL (R) 

SHORT PERIOD (SP) 

SPIRAL(S) 

HEADING (H) 

PHUGOID(P) 

DUTCH ROLL (DR) 

-.849326+00 

-.634506+00 

-.337486+60 

-.2C7046+0C 

-.226786-17 

.101306+00 

•1013OE+0C 

.114516*00 

•114516+00 

I M AG PARTS OF 
0. 

.61 G6NVALUES 
0, 

0. 

0 . 

0. 

.325286+00 

—* 32 5 28 €+00 

.569386+00 

-.569386+00 

EIGENVECTORS 

U cg,h 

V cg,h 

W cg,h 

P h 

q h 

r h 

*h 

®h 


R .300326+00 

.9430«E+0t 

.815236-01 

.683136-01 

—.193066—01 

-.260746-01 

-.910866—01 

,241996-01 

.295646-01 

[-.588936+00 

SP <1 

-.79722E+0C 

-.865126-01 

-.314916—01 

.210296-01 

.378226-01 

*509126—01 

-.360066-01 

-.579396-01 

[-.69560E+00 

-.125826+00 

-.64738 £+00 

.116966-02 

.306746-02 

.010516-01 

. 248376-02 

-.222196-01 

-.269106+00 

S -.269566+00 

, 45238E+OG 

-•69506E— 01 

.630696-02 

-.717326-02 

.171566+00 

-.121236-01 

-. 575366-02 

-.829556+00 

H -.169016-16 

-.125366—16 

-.222016-16 

-.5C138E-19 

• U04 3E-18 

-.226466-17 

.191066-20 

.396106-19 

. 10000E+01 

[ —.673556+00 
[ -.626236+00 

« 11 649E+00 

.298886-01 

.139576-02 

-.106006-01 

-.612406-02 

.135866-01 

-.132536-01 

-.409606-01 

—.369196+00 

• 1C5296— 01 

•4C902E-02 

-.648476-03 

— .1291P C -01 

-.439686-03 

.299526-01 

.443546-02 

| .120766+00 

•40355E+00 

.206636-01 

. 76 0 5 0 E —0 3 

.42716E-02 

-.139036-01 

.589376-01 

-.116656-02 

-.387436-01 

DR<j 









1 • 10955 E— 01 

-.902216+00 

-.194366-01 

.342556-01 

—.431706—03 

-.201576-01 

.997016-02 

-.653866-02 

.169696-01 
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TABLE VI.- Continued 


(b) Values at 30 knots 

TRIMMED VELOCITY - 30.0 KNOTS 7 POINT FORMULA 30 ITERATIONS 


INDEPENDENT VARIABLE TOTAL INCREMENTS 


_l*cg,h_ 
- Tcg,h- _ 

llp.h!l h ' 
'ih 

'2s 

% 

.30488-01 
,30466-01 
. 3048E— 01 
•1745E-02 
.17456-02 
•1745E-02 
.1745E-02 

*IC 

®Oxa 

.17456-02 

.17456-02 

.1745^-02 

.17456-02 

. _ 

.17456-02 



.h 

.17456-02 




A MATRIX 



U cg,h 

v cg,h 

W cg,h 

Ph 

q h 

r h 

» h 

6 h 


*cg,h 

-.28096-01 

-.63598-02 

-.12Z4E-C1 

-.U99E-01 

.84236+00 

.8240E-03 

.26756-07 

-.98026+01 

.90886-19 

ycg.h 

.12656-01 

-.60528-01 

-.90356-02 

-.47136+00 

-.84638+00 

-.15086+02 

.97956+01 

-.10736-01 

.18726—20 

tcg.h 

-.22386+00 

—• 15536-01 

—.55586+00 

—.33346—01 

.15806+02 

-.54096-02 

,37026*00 

* 2840E + 00 

,23246-20 

p h 

.16426-01 

-.15966-01 

-.97016-03 

-.88086+00 

-.13246+01 

.15626+00 

-.16556-06 

•39706-21 

.39706-21 


.17206-02 

•537BE-03 

-.84786-02 

.12226+00 

-.24096+00 

— « 2656E— 02 

•1255E-07 

-.97416-20 

-.97416-20 


-.90556-02 

.23718-01 

•5737E-02 

-.74256-02 

-.43366-01 

-.39496+00 

.35676-06 

-.38126-20 

-.38126-20 

th 

0. 

0. 

0. 

.10006+01 

.10956-02 

-.2897E-Q1 

0. ' 

o. 

0. 

fh 

0. 

0. 

0. 

0. 

.99936+00 

. 3777E— 01 

0. 

0. 

0. 

K 

0. 

0. 

0. 

0. 

-.37796-01 

.99976+00 

0* 

0. 

0. 


B MATRIX 


B IC *TC 9 0m 9 ct ' 

*cg,h .10058+02 *66386-07 -.32566+01 -.11228+00 

*cg,h .13518+00 .98506+01 -.45796+00 ,38198+01 

*eg,h .83646+01 ,30946-05 -.8723E+02 .2839E-01 

Ph .13148-01 . 2 021 E +02 .1162E+01 .17898+01 

4 h -.37508+01 .37476-07 -.93856+00 -.19028-01 

r h -.88496-01 .11936+01 .65106+01 -.36988+01 

$k 0. 0. 0. 0. 

0 . 0 . 0 * 0 * 

'Jijj 0. 0. C. 0. 


REAL PARTS 06 EIGENVALUES 

ROLL(R) SPIRAL (S) HEADING (H) PHUGOID(P) DUTCH ROLL (DR) SHORT PERIOD(SP) 



-.894316+00 

-.512246-01 

,112706-18 

• 62 776 E— 01 

.627766-01 

12986E+00 

-.129866+00 

-.540606+00 

54060 E +00 

IMAG PARTS OF 

0. 

EIGENVALUES 

0. 

0. 

. 286066+00 

-.286066+00 

.668006 *00 

-.668006+00 

.421726+00 

-.421726+00 

EIGENVECTORS 

U cg,h 

V cg»h 

W cg,h 

P h 

q h 

r h 

*h 

9 h 

*h 

R 

-.166316+00 

-.66433E+00 

-.712096+00 

-.955066-01 

.974926-02 

.361256-01 

.107956+00 

-.12419E-01 

-.399716-01 

S 

.718576-01 

-.570186+00 

-• 12758E-01 

.209856-02 

.154056-02 

-.416846-01 

-.645786-01 

.68174E-03 

.814656+00 

H 

.129086-17 

.26637F-17 

-.585476-18 

.326646-2 0 

-.425496-20 

.112576-18 

.187156-16 

,421566-20 

•100006+01 

J 

' -.605106+00 

-.4 11386+00 

.983806-01 

.396946-02 

-.695726-02 

-.740276-02 

-.637146-02 

-. 113476—01 

-.566976-02 

1 

-.565266+00 

-.344476+00 

•12227E+00 

-.283076-02 

-.180886-02 

-.158086-03 

-.160186-01 

.228176-01 

.244076-01 


| .775296-02 

.935236+00 

, 1231 IE— 01 

-.153196-01 

-.865586-03 

.218236-01 

.281396-02 

.172976-02 

-,419516— 01 


[ ,341186-01 

.342416+01 

,497906-01 

-.187026-02 

•21279E-02 

-.247615-01 

,233346-01 

-.275336-03 

-.245546-01 

SP< 

j' -.416446-01 

-.513656+00 

• 284076+C0 

-.241276-01 

.170396-01 

.108986-01 

,390456-01 

-.179906-01 

.76034E— 02 

( -• 40177 E +00 

.816046-01 

-.69323E+00 

.128426-01 

.1483QE-02 

.216785-01 

•78622E-02 

— • 182B9E— 01 

-.34053E-01 
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TABLE VI.- Continued 


(c) Values at 60 knots 


TRIMMED VELOCITY - 60.0 KNOTS 7 POINT FOP »UIA 30 ITERATIONS 


INDEPENDENT variable total increments 


"c€,h • 3O90E-C I 

V cg,K„._ .30986-01 
w cg,h _ .30986-01 

Ph .17956-02 

Jh .17956-02 

T h .17956-02 

• 1795E-Q2 

0^ .17956-02 

.17956-02 


*LC * 17956-02 

A-IC .17956-02 

00m .17956-02 

_ 0 ct _ .17956-02 


A MATRIX 



U cg,h 

V cg,h 

W cg,h 

p h 

q h 

r h 

*h 

9 h 


“cg,h 

-.33306—01 

-.97926-02 

-.19666-01 

— . 859OE-02 

.19396+01 

.20166-02 

.2 71 9 £-07 

-.97956+01 

-.33566-20 

v C g, h 

.69106-02 

-.93906-01 

-.19956-01 

-.16376+01 

-.89696+00 

-.30376+02 

.97926+01 

-.12626-01 

-.90366-21 

Ycg,h 

-.12 5 IE +00 

— • 26966-01 

-.78006+00 

-.52696-01 

.30926+02 

-.11176-02 

•26616+00 

•96966+00 

-.25826-19 

?h 

• 1269E— Q1 

• 9809E-03 

.50336-02 

-.97866+00 

-.13176+01 

.29336+00 

* 1102E-06 

-.79126-21 

-.79126-21 

% 

• 36796-02 

-.19696-02 

-.1098E-01 

• 1196E+00 

-.27236+00 

.13006-03 

♦ 1770E-08 

-.19596-20 

-.19596-20 

Th 

-• 5580E— 02 

.3U9E-01 

.18716-02 

.12886-01 

-.59176-01 

-.55706+00 

-.25616-06 

. 10596—19 

.10596-19 

th 

0. 

0. 

0. 

• 1G0GE+01 

.12896-02 

-.97936-01 

0. 

0. 

0. 

9h 

0. 

0. 

0. 

0. 

.99966+00 

.27166-01 

0. 

0. 

0. 

K 

0. 

0. 

0. 

0. 

-.2719E-01 

.10016+01 

0. 

0. 

0. 





B MATRIX 




B ^C 

6 0m 

9 ct 



u cg ,h 

Y c 8» h 

.10386+02 .50836-05 

-.97176+01 

-.16686+00 



.93706+00 .98006+01 

-.13916+00 

. 9978E+01 



w cg,h 

• 2298E+02 . 1 866E-03 

-.10396+03 

.35796-01 



?h 

-.13226+00 .20176+02 

« 2 8336+01 

.21176+01 



% 

-.35696+01 .25306-05 

-.63936+00 

.27696-01 



I h 

.10226-01 .11926+01 

•5376E+01 

-.95606+01 



th 

0. 0. 

0* 

0. 



Jh 

0. 0. 

0. 

0. 



*h 

0. (. 

0. 

0. 



Rf*AL PARTS 

Of EIGENVALUES 





ROLL(R) 

HEADING (H) SPIRAL (S) 

PHUGOID(P) 

SHORT PERIOD (SP) 

DUTCH ROLL (DR) 


-.956396+00 

-.689866-20 

.565696-01 

.273936-01 

.273936-01 

-.609586+00 

-.609586+00 

-.330156+00 

-.330156+00 

INAG PARTS OF 
0. 

EIGENVALUES 

0. 

0. 

.272896+00 

-.272896+00 

.578616+00 

-.57861E+00 

.988196+00 

—.988196+00 

EIGENVECTORS 

U Cg , h 

V cg,h 

W cg,h 

P h 

q h 

r h 

*h 

9 h 

*h 

R -.799266-01 

-.339286+00 

—.92 5096+00 

-.107356+00 

.359156-02 

.333386-01 

,113906+00 

-.97007E-02 

-.397826-01 

H .222796-18 

-.918806-16 

-• 56071 E-19 

-.326936-21 

.107186-21 

-.688836-20 

-.255356-19 

-.775396-21 

.100006+01 

S .109786+00 

-.72060E+OC 

.289916-02 

-.915896-02 

.100396-02 

-.379016-01 

-•13010E+00 

-.960196-03 

-.670996+00 

( —.926596+00 

pV 

-.292006-02 

.662786-01 

-.129216-02 

— .15230E— 02 

.971296-02 

.85820E-02 

.263396-01 

,999686-02 

[ . 86952E+00 

.217976+00 

• 12992E+Q0 

,259766-02 

.737366-02 

.990956-03 

.623926-02 

.779996-02 

-.169916-01 

f -.293996 +0C 

SP< 

| -.320676-01 

.209996+00 

-.799696+00 

.137966-01 

.299566-02 

.990796-02 

.951196-02 

-.173996-01 

-.17879E-01 

.173896+00 

-.992716+00 

.298506-01 

-.173206-01 

-.123916-01 

-.329316-01 

•125926-01 

•259306-02 

f -.275606-01 

.130596+00 

-.931086-01 

.711376-02 

-.272086-02 

.305206-01 

-.703296-02 

-.673086-03 

— .68989E— 02 

DR< 

[ .619616-02 

•98961E+OG 

.959316-01 

— , 570596-02 

-.199906-02 

.266956-02 

-.338076-02 

.213826-02 

-.286956-01 
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TABLE VI.- Concluded 


(d) Values at 90 knots 

TRIMMED VELOCITY ■ 90,0 KNOTS 7 POINT FORMULA 30 ITERATIONS 


INDEPENDENT VARIABLE TOTAL INCREMENTS 



.30486-01 

?IC . 

,17456-02 

Ycg,h 

.304 8 E-OX 

A xc. 

•17456-02 

*cg,h 

- .30486-01 


,17456-02 

— &— 

.17456—02 

At,__ 

.17456-02 

q h 

.17456-02 




.1745E-02 



—S'- 

_ ,17456-02 




•1745E— 02 



!h_ 

,17456—02 




A MATRIX 



M cg,h 

V cg.h 

W cg,h 

P h 

q h 

r h 


8 h 

♦h 

Vcg,h 

-.42766-01 

-* 5663E— 02 

-.17296-01 

-.94 21 £-02 

.44786+01 

-.14046-04 

,35636-07 

— .97696+01 

-.92946-20 

ycg,h 

■ 5430E-02 

-•1233E +00 

-.22126-01 

-.42386+01 

-.77146+00 

-.45606+02 

. 97646+01 

— .2559E— 01 

-.18596-19 

”cg,h 

— .7645E— 01 

-.41546-01 

-.91476+00 

-.60766-01 

,4606c+02 

.16446-02 

. 29U 6+00 

. 85 82 E +00 

,16526-19 

Ph 

• 1159?— 01 

.15866-01 

.99656-02 

-.96896+30 

— .12456+01 

,30708+00 

,29396-06 

,29656—20 

.29656-20 

$h 

• 302 5 r -02 

-.3043.--02 

-♦U29E-G1 

, 1174E + O0 

-.26246+00 

.73566-0? 

-.4O26E-08 

-• 6353t-2i 

-.63536-21 


—,4397c -02 

.3452f-»U 

. 53266-02 

, 2505 E— 01 

-.52396-0! 

-.65056+00 

-.t670E-0fc 

-* 18631-19 

-.13636- IV 

*h 

0, 

n. 

C. 

.10006+01 

,26196-02 

-.5785E-01 

0. 

0. 

0. 

9h 

c. 

1 * 


U. 

.99961+00 

.29806-01 

0, 

0. 

0. 

K 

0. 

0. 

0. 

0* 

-.29926-01 

.10036+01 

0. 

0, 

u. 



ROLL(R) HEADING (H) SPIRAL (S) PHUGOID(P) SHORT PERIOD (SP) DUTCH ROLL (DR) 
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Figure 1.- Simulation mathematical model block diagram. 









Figure 1.- Continued. 
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Figure 1.- Continued, 
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Figure 1.- Continued. 
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Figure 2.- U.S. Army CH-54 helicopter and concrete block load, showing helicopter 
body axes, main- and tail-rotor shaft axes, and Earth axes. 






Figure 4.- Load suspension cable angle definitions. 
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Figure 12.- Simulation control console 
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Figure 16.- View of simulator cockpit showing cyclic, collective, 

and pedal controls. 
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Figure 17*- Simulator cockpit showing pilot's visual display monitor 
lens systems and researcher's black and white monitor. 





Figure 18. - Simulator RESET, HOLD, and OPERATE push buttons as well as 
AFCS, YAW, altitude (ALT), and HOVER push buttons and indicators. 
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Figure 19.- Simulator control system analog computer 

(on the left) . 
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(a) General flow diagram. 


Figure 21.- Simulation flow diagram 
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(b) Mathematical model flow diagram. 
Figure 21. - Concluded. 
















L-75-7494 


Figure 22*- VLDS model board mounted vertically 
and camera drive system* 
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Figure 23.- VLDS heliport landing zone close-up 
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Figure 24. ~ Sketch of computer-generated load/landing 

zone display. 
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Figure 26.- Photograph of pilot's out-the-window display at an 
altitude of 55 m and a distance of 30 m from landing zone. 





Figure 27. ~ Photograph of pilot’s out-the-^ 


\o 

CM 

CM 

in 

r 

r** 

r- 

r 

a 


£ 

rtf 

£ 

V0 

vO 

H-J 

0 

a) 

£ 

-P 

*H 

■P 

r— I 

rtf 



Time, sec 


Figure 28.- Comparison of flight data with simulation data 
at hover for identical control inputs. 
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Figure 29.- Comparison of flight data with simulation data 
at 38 knots for identical control inputs. 



rad 0 u ,rad 






1. Report No, 

NASA TP- 12 85 

2. Government Accession No, 

4. Title and Subtitle 

DEVELOPMENT AND VALIDATION 
A HELICOPTER AND EXTERNAL , 

OF A PILOTED SIMULATION OF 
SLING LOAD 

7. Author(s) 

J. D. Shaughnessy, Thomas 
and Kenneth R. Yenni 

N. Deaux, 


9. Performing Organization Name and Address 

RASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 


National Aeronautics and Space Administration 
Washington, DC 20546 


3. Recipient's Catalog No. 


5. Report Date 

January 1979 


6. Performing Organization Code 


8. Performing Organization Report No. 
L- 11925 


10. Work Unit No. 

505-08-26-01 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 


15. Supplementary Notes 

J. D. Shaughnessy and Kenneth R. Yenni: Langley Research Center. 

Thomas N. Deaux: Sperry Support Services, Hampton, Virginia. 


16. Abstract 


A generalized, real-time, piloted, visual simulation of a single-rotor helicopter, 
suspension system, and external load is described and validated for the full 
flight envelope of the U.S. Army CH-54 helicopter and cargo container as an 
example. The mathematical model described uses modified nonlinear classical 
rotor theory for both the main rotor and tail rotor, nonlinear fuselage aerody- 
namics, an elastic suspension system, nonlinear load aerodynamics, and a load- 
ground contact model. The implementation of the mathematical model on a large 
digital computing system is described, and validation of the simulation is dis- 
cussed. The mathematical model is validated by comparing measured flight data 
with simulated data, by comparing linearized system matrices, eigenvalues, and 
eigenvectors with manufacturers 1 data, and by the subjective comparison of 
handling characteristics by experienced pilots. A visual landing display system 
for use in simulation which generates the pilot* s forward-looking real-world dis- 
play is discussed, and a special head-up, down-looking load/landing zone display 
is described. 


17. Key Words (Suggested by Author(s)) 

Helicopter 
External/ sling load 
Rotorcraft 
Visual simulation 


19. Security Classif. (of this report) 
Unclassified 


18. Distribution Statement 
Unclassified 


Unlimited 


20. Security Classif. (of this page) 
Unclassified 


Subject Category 05 


21. No. of Pages 
111 


22. Price* 
$6.50 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1979 















National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business 

Penalty for Private Use, $300 


THIRD-CLASS BULK RATE 


National Aeronautics and 
Space ^Administration 







v- : : . . ; ; '. . gf| 


: < ' : 

1 1 i ^ * V i ^ liS H ,* s i >« f SIB H lilp M i irtilil MMMiife 






POSTM ASTER * If Undeliverable (Section 1 58 

Postal Manual) Do Not Return 



